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Uber die relativistische Elektronenoptik elektrostatis cher 
: Beschleuniger 


Von Ents B. Bas, LucreN Preuss und WERNER SCHNEIDER, Ziirich?) 


1. Einleitung 


Ein elektrostatischer Partikelbeschleuniger besteht aus einer Kette Pie 
Beschleunigungsstufen, zum Beispiel in der Form von Lochblenden- oder Rohr- 
_ Immersionslinsen (Figur 1). Die Aufgabe dieser Linsenkette ist erstens die 


wi ae; He 3 ty ty kg kp Hg xy Ky Zq-6050m 
450m 
, 405¢m —=— 200¢mM : 
L. 605¢em = 
Figur 1 


Prinzipschema eines elektrostatischen Beschleunigers (Abmessungen gelten fiir ein praktisches 
Beispiel). x, = U,/U, usw. 


-Beschleunigung der von der Partikelquelle A emittierten Ionen bzw. Elek-. 

' tronen auf die Endenergie Up. Zweitens muss sic den schmalsten Querschnitt 
des Ionen- bzw. Elektronenstrahles am Eingang im Punkte P auf das Target 
im Punkte Q abbilden. Die Elektronenoptik dieser Abbildung ist der Inhalt 

dieser Mitteilung. Die Aufgabe ist, dabei die Beziehung zwischen den Achs- 

_koordinaten der beiden konjugierten Punkte P und Q aufzustellen, wobei 

noch die laterale Vergrésserung f und die Winkelvergrésserung y anzugeben ist. 
Fiir die achsennahen Elektronen- bzw. Ionenbahnen in einem rein elek- 
trischen Feld gilt allgemein folgende Differentialgleichung®) : 


ar ae ae 


/(L+e®)o [Va+ea)® =| ee ; fl 1 22O)O"7=0. (1) 


1) Institut fiir technische Physik der Eidgendssischen Technischen Hochschule. 
2) W. Guaser, Grundlagen der Elektronenoptik (Springer, Wien 1952), S. 144. 
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@ bezeichnet das elektrische Potential entlang der Achse, und (1 + e@®) 


baw. (1 + 2 ¢®) ist das relativistische Korrekturglied mit 


pS oo) ele (2) | 


2 Mg 6 


e Ladung des Elektrons bzw. Ions; 
my Ruhemasse des Elektrons bzw. Ions; 
c  Lichtgeschwindigkeit. 


Wir erhalten 
fiir Elektronen ¢ = 0,978-10-§ [V4]; 


fiir Protonen e = 0,539-10-® [V4]. 


Fiir das relativistische Korrekturglied 1 + ¢ ® erhalt man somit zum Beispiel | 


bei einem Potential von ® = 5-108 V 
fiir Elektron 1 + ¢® = 5,88; 
fur Proton © 1 +e@ = 10027. 


Da elektrostatische Beschleuniger nur bis zu einer Spannung von 
U,, x 5-108 V gebaut werden, kann bei den Ionenbeschleunigern die Relativitat 
stets unberiicksichtigt bleiben. Hingegen sehen wir, dass. bei den Elektronen- 
beschleunigern der Relativitatsfaktor betrachtliche Werte annehmen kann, so 
dass es sich hier aufdrangt, itber dessen Einfluss Rechenschaft zu geben. 


2. Relativistische Linsendaten einer elektrischen 
Zweirohr-Immersionslinse 


Ein elektrostatischer Beschleuniger mit niedriger Stufenzahl, wie es fiir die 


Freiluftausfiihrungen tiblich ist, besteht aus einer Kette von sogenannten — 


elektrischen Zweirohr-Immersionslinsen. Eine typische Geometrie einer sol- 


chen Zweirohrlinse ist in Figur 2 dargestellt. Die elektrischen und geometri- 


schen Parameter der Linse findet man in Figur 2 ebenfalls angegeben. Die 
Kardinalelemente der Linse werden in dieser Arbeit wie folgt bezeichnet: 


Bildseitig: Brennweite = /’; 
auf die Mittelebene bezogene Brennpunktkoordinate zp, = F’; 
auf die Mittelebene bezogene Hauptebenenkoordinate z,,, = H’. 


Dingseitig: Wie oben, jedoch ohne Strich. 
U’ bezeichnet die bildseitige und U die dingseitige Elektrodenspannung. 


x = U/U’' ist der elektrische Linsenparameter; da es sich um eine Beschleuni- 
gungslinse handelt, ist stets 


Urea aeallk | 


oa ‘a . eos . ~ y 
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Die Brechkraft einer nicht relativistischen’ Linse wird durch die beiden 
Parameter, den elektrischen Parameter und dén geometrischen Parameter 
S/R, bestimmt. Bei den relativistischen Linsen kommt noch die Abhangigkeit 
von U' hinzu. In dieser Arbeit werden Wir uns tibersichtlichkeitshalber nur auf 
einen Wert des Geometrieparameters, 

a 


R = 08 , 


beschranken. Die Ergebnisse fiir andere S/R-Werte werden in einer spateren 
Arbeit publiziert. 


Figur 2 
; Geomeirie einer typischen elektrischen Zweirohr-Immersionslinse fiir elektrostatische Beschleuniger. 
Brennweiten: f = 2(H) — 2(F); ff = 2(H’) — 2(F’); f = Vets iO el 


Zur Berechnung der Kardinalelemente der Linse muss man als erstes den 

Verlauf des Achsenpotentials ®(z) kennen. Wir haben diesen Verlauf in einem 

- Widerstandsnetzwerk-Modell nach LrzBMANN3) ermittelt. Die Maschenweite 

im Netzwerkmodell entsprach 1/20 R. Der Messbereich entlang der Linsen- 

_achse betrug 160 Maschenweiten. Die Potentialverteilung entlang der Linsen- 

_achse wurde nur einmal fiir die Elektrodenspannungen 0 und 1 V gemessen. ’ 
- Bezeichnen wir dieses Potential mit w(z), 


0< plz) <1, (3) 


so ergibt sich bei einer vorgegebenen x und U’ fiir den tatsdchlichen Potential- 
_verlauf auf der Linsenachse 


: (2) = (1—) U' [-*— + ve]. (4) 


Die Berechnung von @(z) nach dieser Beziehung wurde in das — im folgenden 
noch zu besprechende — Rechenprogramm einer elektronischen Rechenmaschine 
 eingebaut. 

< 3) G. Liepmann, Solution of Partial Differential Equation with a Resistance Network Analogue, 
; Brit. J. appl. Physics 1, 92 (1950). 

A 
4 
“ 


fo 
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Zur Berechnung der Kardinalelemente einer elektrischen Immersionslinse 
miissen zwei Elektronenbahnen im Linsenfeld durchgerechnet werden, das 
heisst, wir miissen zwei partikulare Lésungen der Differentialgleichung (1) 
bestimmen. Hierfiir ist es am zweckmassigsten, zwei besondere Elektronen- 
bahnen zu betrachten: eine von links parallel zur Achse im Abstand 1 ein- 
fallende Bahn o0(z) und eine von rechts parallel zur Achse im Abstand 1 ein- 
fallende Bahn o(z) (siche Figur 3). Bezeichnen wir die Asymptoten dieser 


| 9a 


Figur 3 


Hilfsbahnen fiir die numerische Berechnung der Kardinalelemente einer Immersionslinse. 


Bahnen nach dem Verlassen des Linsenfeldes mit g(z) bzw. h(z), so sind die 
asymptotischen Kardinalelemente der Linse gegeben durch: 


Fi = %, (5) 
f = sar 6) 
ory a 
i= anaz a 


7 
wobei 22 und z} die Schnittpunkte der Geraden g bzw. h mit der Z-Achse 
bedeuten. | 

Die Differentialgleichung der Bahnen in der Form von Gleichung (1) ist — 
fiir die numerische Integration nicht sehr geeignet, da dort die zweite Ans 
leitung von (z) auftritt, was eine zweimalige numerische Differentiation von 
(z) verlangt. Es ist deswegen zweckmissig, nach Picut die Gleichung (1) wie 
folgt umzuformen: Man fiihrt eine neue Bahnvariable R(z) ein, nach der 
folgenden Transformationsgleichung: 


¥(2) = 7 ————__— R(z). (9) 
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Fiir R(z) gilt dann folgende Differentialgleichung: 


a eae a: : 
RG eg) R=0 (10) 
mit | 
| 1+ = (1+) 
AO) erat (11) 


Bei den nichtrelativistischen Linsen ist G = 1, und die Bahngleichung erhalt 
folgende Form: 


1 
ney = v5 Re) (12) 
‘ a @’ \2 
Rey re le) R= 0. (13) 


Fiir das Verhaltnis zwischen der relativistisch gerechneten Bahn rp(z) und 
der nichtrelativistisch gerechneten 7y(z) kénnen wir damit schreiben: 


re(2) 1 Rez) _ oe Ree) 
rx@) — Yire® Rye) ~ * Rey ee 
mit 
1 
ASS oD) es A ee 


Rp(z)/Ry(z) wird nur durch den Verlauf von G innerhalb der Linse beein- 
flusst, das heisst, das Verhaltnis zwischen relativistischen und nichtrelativi- 
stischen Bahnen wird durch Einwirken zweier Faktoren K und G (besser 
gesagt durch 1/G) bestimmt. Wir werden darauf noch zuritickkommen. 

Die Integration der Bahngleichung (1) bzw. (10) wurde nach dem Ver- 
fahren von RunGE-Kutta mit Hilfe einer elektronischen Rechenmaschine 
durchgefiihrt. Die ®-Werte wurden aus den im Netzwerk gemessenen Werten 
nach Gleichung (4) von der Rechenmaschine ermittelt. Dabei wurden nach der 
kubischen Interpolationsformel noch Zwischenwerte berechnet, so dass fiir die 
Integration der Rechenschritt auf die halbe Maschenweite (1/40 des Rohr- 
radius der Linse) reduziert wurde. Das automatische Rechenprogramm der 
elektronischen Rechenmaschine umfasst die Berechnung von @(z) nach 
Gleichung (4), Interpolation, Integration der Gleichung (10), Riicktransforma- 
tion nach Gleichung (9) und die Berechnung der Kardinalelemente der Linse 
nach den Gleichungen (5) bis (8), so dass die Kardinalelemente direkt aus den 
- Registrierstreifen der Maschine abgelesen werden konnten. 

In Figur 4 ist das Verhaltnis zwischen relativistischen und nichtrelativisti- 
schen Koordinaten der bildseitigen Brennpunkte Fp/fy in Abhangigkeit von 
U’ mit x als Parameter eingetragen. Dort findet man auch’ versuchsweise 
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K(®) und 1/G(®) fiir das bildseitige Potential U’ berechnet und eingetragen. 
Der Einfluss von K und G auf den Verlauf der Kurven ist deutlich. Wahrend 
bei den kleineren x-Werten, das heisst starkeren Linsen, hauptsachlich G den 
Kurvenlauf bestimmt, zeigt sich bei schwacheren Linsen immer mehr der 
Einfluss von K. In Figur 5 findet man die entsprechende Darstellung fiir 


i gen T T T tT — T 1 re ae | 
500 7000 1500 2000 

U'(kV] —— 

Figur 4 


Bezogene bildseitige Brennpunktlage der relativistischen Zweirohr-Immersionslinse in Abhangigkeit 
vom bildseitigen Linsenpotential U’. 


frify. Auch diese Kurven haben ganz ahnlichen Verlauf wie die Kurven fiir 
Fp/Fy. Auffallend ist noch, dass die Kurven fiir sehr starke Linsen (% = 0,1 
und 0,2) durch die Uberschneidungen einen riicklaufigen Charakter aufweisen. 
Dies ist in der Darstellung in Figur 6, wo f/f in Abhangigkeit von x bei 
konstantem U’ aufgetragen ist, besser zu itibersehen. Bei kleineren Linsen- 
spannungen hat x sehr geringen Einfluss auf f/f. 

Fir die praktische Beniitzung der wiedergegebenen Rechnungsergebnisse 
ist es meistens giinstiger, anstelle von Brennpunktkoordinaten von den Haupt- 
ebenenkoordinaten H und H’ auszugehen; es gilt dabei 


H' — a le jae (16) 


Hees | eee Jets 


rer). ee 


~*~ 


Pe ee ee ee ee Se et RT 


r ¢ ‘aa = . tea’ a4 ~ 
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FO id 71000 7500 = =——«000 
Figur 5 
Bezogene bildseitige Brennweite der relativistischen Zweirohr-Immersionslinse in Abhangigkeit 
vom bildseitigen Linsenpotential U’. 


0 a5 10 


Figur 6 


Bezogene bildseitige Brennweite der relativistischen Zweirohr-Immersionslinse in Abhangigkeit vom 
Spannungsverhiltnis ¥ = U/U’. 
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In Figur 7 findet man Hp — Hy und in Figur 8 Hp — Hy in Funktion 
von U’ mit x als Parameter graphisch dargestellt. Eine Wiedergabe von fz//y 
eriibrigt sich, da fp nach der Beziehung: | 


U (1+ €U) , 
Ir =F (1+ ¢U’) Tr ue 


aus f berechnet werden kann. Damit sind wir in der Lage, aus den nicht- 
relativistischen Linsendaten die relativistischen Linsendaten der Rohrlinse 
S/R = 0,8 fiir alle praktisch interessierenden Werte von x und U’ zu ermitteln. 
Die Zahlenwerte der Kardinalelemente /, /’, H und H’ einer nichtrelativisti- 
schen Zweirohr-Immersionslinse mit S/R = 0,8 sind in der Tabelle 1 zusammen- 
gestellt. Die Elektronenbahnen in den nichtrelativistischen Linsen sind vom 
Verhaltnis der angelegten Elektrodenspannungen und nicht von ihren Absolut- 
werten abhangig, so dass in der Tabelle 1 nur x als elektrischer Parameter 
auftritt. Sowohl die Daten in der Tabelle 1 wie die Angaben in den Figuren 4 
bis 8 beziehen sich auf eine Beschleunigungslinse. Sie konnen selbstverstand- 
lich auch fiir eine Verzégerungslinse beniitzt werden, wenn sinngemass die 
Bezeichnung «bildseitig» und «dingseitig» vertauscht wird. 


Tabelle 1 


Kardinalelemente dev nichtrelitivistischen Zweivohy-Immersionslinse mit S/R = 0,8 
(4 = U|U’, Beschleunigungslinse) 


x | 0,0666] 0,1 0,2 | 0,333 | 0,5 0,666 | 0,75 0,85 0,95 0,9666 
i 
R 1,312} 1,865] 4,103] 9,468] 25,40 Visas)  L55,61) 48a, 8 Loy 8473 
h 
RP —4,818 | -5,687 | —9,009|-16,23 |-35,71 |-94,64 [|-179,3 |-521,9 |-1166 8616 
H 
R 1,547 |-1,651 | -2,060|— 2,832|- 4,387|- 7,447 |- 10,42|— 17,93|— 26,52|- 64,78 
H’ 
= —2,680 | —2,712 |-3,022|- 3,740 |— 5,261|- 8,295 Li25 |= Uss7Z 27,27\— 65,37) 


3. Relativistische Elektronenbahnen 
in einem praktischen Beschleuniger 


Wie im vorigen Abschnitt gezeigt wurde, kénnen sich die relativistischen 
und die nichtrelativistischen Brennweiten bis zu 20°, voneinander unter- 
scheiden und die relativistischen Hauptebenen bis zum fiinfzehnfachen Wert 
des Rohrradius von den nichtrelativistischen Hauptebenen entfernt liegen. Es 
ist nun interessant zu untersuchen, wie sich diese Differenzen bei der Berech- 


= 


es | 


* 


—_——— —_=" = =— =) on 


 . 
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Figur 7 


Bezogene Hauptebenenlage der relativistischen Zweirohr-Immersionslinse in Abhangigkeit vom 
bildseitigen Linsenpotential U’. 


ee a = 
0 500 7000 1600 2000 
UK) —— 
Figur 8 


Bezogene dingseitige Hauptebenenlage der relativistischen Zweirohr-Immersionslinse in Abhangig- 
keit vom bildseitigen Linsenpotential U’. 
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nung eines mehrstufigen Elektronenbeschleunigers auswirken. Wir betrachten 
dafiir einen praktisch ausgefiihrten 10stufigen Elektronenbeschleuniger fiir 
1500 kV. Die Linsenlagen und die Linsenspannungen bzw. x-Werte sind in 
Figur 1 angegeben. Das Target Q befindet sich 650 cm von der als Koordinaten- 
ursprung gewahlten Mittelebene der ersten Linse entfernt, also 


to = 605 cm. 


Der Radius R der Rohrlinsen betragt 5 cm. Der x-Wert der ersten Linse , ist 
als verdinderlich anzusehen, da die Wahl der Elektroneneinschussenergie, das 
heisst der Anodenspannung U, der Elektronenkanone, noch freisteht. Wir 
werden also die Lagekoordinate z, des Dingpunktes P sowie die Lateral- 
vergrésserung f und die Angularvergrésserung y in Abhangigkeit von x, (bzw. 
der Anodenspannung U,) zu berechnen haben, wobei noch zu beachten ist, 
dass zwischen f und y folgende Beziehung besteht: 


i Uy dere 

sosiok) \ Un i re 0) ; — 
In Figur 9 ist das Ergebnis der relativistischen und nichtrelativistischen 
Rechnung graphisch dargestellt. Wie man sieht, ist der relativistische Effekt 
im Endresultat nicht sehr gross. Der maximale Fehler in der Angabe von 2, 
betragt nur 11 cm bei x, = 1. Auch die Lateralvergrésserung wird gegeniiber 
dem relativistischen Fall maximal um etwa 40° zu gross berechnet. Interes- 
sant ist es, dass durch die Beziehung (19) diese Differenz der Lateralvergrésse- 
rungen wieder kompensiert wird, womit der Unterschied bei den Angular- 

vergrésserungen fast vollstandig verschwindet. 
In Figur 10 ist fir Uy =50kV die Elektronenbahn durch den ganzen 


Beschleuniger dargestellt. Die Bahn ist riickwartsgerechnet. Der Winkel gegen 
die Achse am Ausgang des Beschleunigers betrug in beiden Fallen 


05 = IO . 


Zusammenfassend kann gesagt werden, dass in den meisten Fallen die 
Elektronenoptik eines elektrostatischen Elektronenbeschleunigers nichtreiativi- 
stisch berechnet werden kann. 

Fiir die numerische Rechnung wurde die elektronische Rechenmaschine — 
ERMETH des Institutes fiir angewandte Mathematik der Eidgendéssischen 
Technischen Hochschule verwendet. Wir méchten auch an dieser Stelle dem 
Leiter des Institutes fiir angewandte Mathematik, Herrn Prof. Dr. E. STIEFEL, 
fiir sein freundliches Entgegenkommen sowie Herrn Prof, Dr. H. RUTIS- 
HAUSER und seinen Mitarbeitern fiir ihre wertvollen Ratschlage unsern besten 
Dank aussprechen. Wir danken auch Herrn Prof. E. BAUMANN, Leiter des 
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Institutes fiir technische Physik der Eidgenéssischen Technischen Hochschule, 
fiir die Erlaubnis zur Veréffentlichung dieser Arbeit. 
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Figur 9 
Ergebnisse der relativistischen und nichtrelativistischen Berechnung eines elektrostatischen 10stu- 
figen Elektronenbeschleunigers. z,, Dinglage bei einer Bildlage z9 = 605 cm (bezogen auf die Mittel- 
ebene der ersten Beschleunigungslinse), / Lateralvergrésserung, y Angularvergriésserung, Uy Ein- 


schussenergie der Elektronen. 1500 kV, S/R = 0,8, R = 5cm. relativistisch; ——— nichtrelati- 
vistisch. 
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Figur 10 


Rine relativistisch und nichtrelativistisch berechnete Elektronenbahn in einem 10stufigen Elek- 
tronenbeschleuniger. Uy = 50 kV; Ug = 1500 kV; —— relativistisch,; ——— nichtrelativistisch. 
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Summary 


Relativistic and nonrelativistic data for the coaxial two-cylinder electrostatic 
electron-lens are compared. Values of the potential along the axis were measured on 
a resistance network analogue, whilst trajectories were calculated on a high speed 


digital computer. pine 
Index FR indicating relativistic values, and N non-relativistic ones, quotients 


Fi/Fy of the focal point coordinates (image sided), frit of the focal length 
(image sided) and differences of the principal plane coordinates (HR — Hy)/R (ob- 
ject sided) and (Hz — Hy)/R (image sided) are represented graphically as functions 
of the imagespace potential and U’. The gap between the cylinders was In all a 
0,8 R (R radius of the cylinders). U being the potential on the object side, and U on 
the image side, x = U/U’ was varied from 0,1 to 0,966 and U’ from 0 to 2,000 kV; in 
the relativistic case a maximal departure of about 20% from the nonrelativistic 
values was found for focal length. A ten-lens accelerator of 1,5 MeV overall tension 
has been calculated both ways. 


(Eingegangen: 16. Mai 1959.) 


Reverse iv situ Combustion 
By Rosert P. Girpert, Pittsburgh, Pa., USA) 


Introduction 


In this paper we are going to investigate the 7 sitw combustion of crude 
oul, that is the ‘burning in place’ of oil which has saturated a porous rock 
formation. This phenomenon is of interest to the petroleum industry as a 
thermal-recovery method for residual oil in a reservoir; the underlying concept 
of all such methods is that by raising the temperature of the reservoir, the 
viscosity of the oil is lowered and thereby its mobility is increased. 

Usually, 2m situ combustion is instigated by first heating a local region and 
then supplying sufficient oxygen to support combustion. The physical phe- 


nomena which then take place are quite involved, for instance, heat is transferred 
by three different mechanisms: 


(1) Convection by gases which flow through the porous rock carrying their | 


heat content with them, 

(2) conduction in the oil-bearing rock formation, and 

(3) exchange of heat between the rock formation and the gas. 
Hydrodynamic effects, such as varying oil-water-gas ratios, condensation of 
hot vapors outside of the combustion zone, and the gravity segregation of 
various fluids, are usually influential factors in burning. Furthermore, the 


1) Mathematics Department, The University of Pittsburgh. 
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heat-source describing combustion is a non-linear function of the temperature 
and the concentration of the reactants. 


In order to simplify the problem somewhat, we are going to investigate the 
case of tar-sand reservoirs. Here, thé initial viscosity of the oil before heating 
is so high that it may be considered immobile. 


Mathematical Formulation of the Problem 


It would be well to list here the symbols used in the following work. 


Variables 
+ distance along the direction of burning; 
t time; 
U, combustion zone velocity; 
Vv, gas velocity; 
2=x-—v,t distance measured in a coordinate frame moving with the combus- 
tion zone; 
? temperature; 
F. temperature of the solid matrix; 
i temperature of the gases; 


p partial pressure of oxygen; 
Rip, T) an arbitrary heat-source function, measured in units of mass hydro- 
carbon consumed per unit of time. 


Constants 
f porosity of the solid; 
0, density of the solid; 
0, density of the gases; 
00 density of the oil; 
e, heat capacity of the solid; 
Cy heat capacity of the gas; 
oy heat capacity of the oil; 
ee residual oil saturation ; 
P. partial pressure of the oxygen before entering the combustion zone ; 
me ATT? the heat of combustion of the oil in units of heat per unit of mass; 
kK conductivity of the solid matrix containing the oil and gas; 
Bk diffusivity of the solid matrix; 
h, heat exchange coefficient between solids and gases; 
LE the steady-state combustion zone length; 
1 the maximum temperature ; 


a; the ambient temperature. 
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Let us consider the situation where oxygen is supplied to the reaction zone 
from its cold side; in this instance because of the immobility of the cold oil, 
any hydrodynamic effects must occur in the combustion region itself, or in that 
region through which it has just passed. We shall refer to this situation, where 
oxidant and reaction zone proceed in opposite directions, as reverse im situ 
combustion. ax 4 

Now, let us investigate reverse im situ combustion occuring in a one-dimen- 
sional model of a porous medium, where the pores are partially filled with an 
immobile hydrocarbon. Let us assume that the rate at which hydrocarbon is 
consumed at any one point is described completely by a continuous function 
of p and T, where # and T are measured at that point. Furthermore, let us 
assume that the important mechanisms of heat transfer involved here are (1) 
convection and (2) conduction. In this case, within the combustion zone, the 
differential equation, 


{(1 — f) 0,6 + FIL — $9) ¢, @ + So Go Col} = 1,0, es (c, 0, T) | 
gar (2 ee 


governs the transfer of heat. If the heat content of the oxidant and that of the 
products of combustion are relatively the same, then 


0 BIE 
Ox (c, Qs Le hed Og Os oy 
and 
on Oe BORE 
sealers fly Og Vg cate LOR by 1.) ae roe (2) 
where 


a= {(1 — f) O;, 9  aaSq) Op fz + Sq Qo Col} - 


We introduce the heat diffusibility of the solid matrix as K /« = k measured 
in length? per unit time, and the ‘velocity’ w = f Cg Og V,/%; equation (2) then 
becomes 

GUE 

oye (2a) 


By assuming a steadily-propagated combustion exists, that is a steady-state 
condition in a coordinate frame moving along with the combustion zone, we 
can simplify this problem considerably. If distance in the combustion zone is 
measured in terms of z, where z = x — v, t, then equation (2a) can be written as 


Oh oT AH® 
i, oe 


aa LE ey @T 
OT Be Saar Teg LO lis raat (3) 


If the rate at which heat is produced is proportional to the rate at which 
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oxygen 1s consumed, then for any point in the reaction zone we may write 


q - 
<P 0, Rip, T), (4) 


where Cy is a constant of proportionality. We can determine this constant, by 
integrating equation (4) over the combustion-zone and introducing the bound- 
ary condition p = P, when z = L. It follows that 


“ 


and that 


Now, since 


hen | 
al 


Rh 
= 
> 
3 

> 


se Sal 


is the total rate of oxidant consumption in the reaction zone, this quantity 


! must be equal to the influx of oxidant in the moving zone, that is 

4 

¥ L 

° ’ 5 

, f (0, +2.) = | RUB, T) de. 

i 0 

_ We obtain from this the following integral relationship 

P, P, , | 

See LA Pee 5 
PP 7G a fs (6, 7) (5) 


Our problem now, is to solve equations (3) and (5) simultaneously. By intro- 
ducing a change of dependent and independent variables we are able to reduce 
‘these equations to a single first-order differential equation. We proceed as 


follows by letting 


(ie) *=a, 6 


whereby 
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(It is assumed of course that all necessary terms are monotonic.) Equation (3) 


becomes 


AH? dp dp 
(Ot ONO a a 
or 
AH® 
—(w+v,) aT + —— t, dp=hkdg, (7) 


which integrates into 


AH® 
— (wat) Te Cp as ae (8) 


c, is a constant of integration, which may be determined from a boundary- 
condition such as T = J, and wy = 0 at z = L; this leads to an equation of the 


form 
AH 


lee 


w+) (T-T) +7 2,6 -p) = he. (9) 


Making the previous change of dependency in the differential form of equation 
(5), and substituting the result of equation (9) into this yields 


AH® 
aD. G06 Pe 
dp Ris Take R@T) 


An equation of this form frequently may be solved by an iteration method; 


the requirements necessary for a solution are known as the Cauchy-Lipschitz 


conditions. 
Consider the general, first-order, ordinary, differential equation 
or in 


op = Fe. T). 


A unique, continuous, solution exists for this differential equation, which ~ 


passes through the point (f,, 7;), if there is a region containing this point in 
which F(p, T) is continuous, and if the following requirements hold: 

(Ly) JP (Dp . Ti) = BD. Ts) | Sb T, Sl ewhere tf asatine Lipschitz constant. 
(2) |F(p, T)| SM, that is F is bounded in this region. 


(3) We can bound our region by the inequalities \p—,| SA and © 


fies. 


Referring back to equation (10), we notice that if R(p, T) + 0 anywhere 
in the combustion zone our requirements are satisfied. We now define the 
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functions 
p 7 
T,(6) = T+ | Fb, Ty) dp, 


v 
?, 0 


p 
T,(d) = Ty+ | Fb, Ti) 4p, (11) 
Po 
p 


T,(0) = Ty + | F(p. T,_.(0)) dp . 


Po 


By the Cauchy-Lipschitz theorem these successive approximations to the solu- 
tion of equation (11) converge, and furthermore the function 


T'(p) = lim T,(6) 


n—>co 


is actually the solution to our differential equation. 
Now, assuming that a solution has been found we have 


dp adi \-1 dT 
“dz = (ip) dz = ¢, R(f, IVs 
so that one may write 
: or dT 
r =i 
=a) (ae) Rp. Ton ve 


Equation (12) is the steadily-propagated temperature profile. The cor- 
responding oxidant partial-pressure profile may be obtained from 


P 
dp 
= “| R(p. T(P)) * ee 


In all of the indicated computations, one must be aware which functions 
are monotonic, and be careful to choose the correct branches of non-monotonic 
functions. 

Next, we turn to the situation where there are three modes of heat transfer 
(1) convection, (2) conduction, and (3) exchange between gases and solids. 
The differential equations which govern the transfer of heat in the reaction 
zone for the gaseous and solid media are respectively 


oT, ee | Y - 
FL 0, GE = Rb, T) — bE, — TF) (14) 


_ -ZAMP X/35 
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and 


Or age 
8 Sia a 1 
at k axe h(T, T,) . ( 5) 


Equations (14) and (15) obey a general equation of the form 


iy SE +o aa Bea ena ee F: 
A soe cy) 5) = y(p, T), 
| Re, Lj wren. JL. 
where y(p, T) = ) a, Rie r) + s i: nee Ce ee 


As we did earlier, we again transfer coordinates to the moving frame z = x — 9, f, 
and make use of our change of dependency 


-( 0% \-1 
a 


> (4) = (7), 


eT a2T Ze 

Ox? )> de ere 

as aT Be dee 
— 2a saa me le 

mae) as. © ae + Pl) , 


 Oaf0T ee Oo 


( 
( 
ie a _ __(ex/o2) 
| 


OFT dy \3 
oe Berry la ; 


Upon using these transformations and the notation yp =dg/aT, equation (16) 
becomes 


% (he +h) p— uv py +k lv. g (pe) +», oe’) — PY] 


; (17) 
+ U, Lu, p (p ¢’)’ — h, |] = w(p, T). 


For the case T= 2 w(p,. 2). h, R(p, T,), and as before we have 


dp dp 
Gie tedr ea R(p, ib: 
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Equation (17) then becomes 


U, (h, oe h,) re v, re +k [(u, te V,) Ys ®, = h, P<] = Ug U.(P, p.)’ i Ug h, = Sy es (18) 


which integrates into 


[v, (/, a h,) ae, Ue h, ge cn ve Ps ae k [(v, is 2,) Ps , = he Ps! | 
; feasted 
a U0, PO, = C, Dar Ca, | 
or 
aes _ — [v, (Ag + hg) — he) 7, + (5 = & hh,) + ¢,% oe 
ebb Pvp.) =.= (hs + hg Ges es c g) Ps iP + & 
Ps (p s Ps) [Rk (uv, + Vg) aD, V,) Ps : 
However, from above we have ¢, ~’ = c, R, so 
, , i f OR oR , 
P,P + P,P Sere ot Op ') =, Rp, DON 
or 
P' F(p, TG.) + 9B = 4 AE, T) - (19) 


We may now remove @,, by substituting in 


R(p, Ts) 


2 a po ; 


from this we are led to an equation of the form 
Pe GP pol), 


which may be solved by reducing it to two first-order simultaneous differential 
equations 
Piss beet Shia eT 


Zusammenfassung 


Fiir den Fall einer hinsichtlich Temperatur und Partialdruck des Oxydenden 
nichtlinearen Quellenfunktion wird die in-situ-Verbrennung von Rohdl untersucht. 
Die Warmeiibertragung wird in ihrer Abhangigkeit, von den Mechanismen der Kon- 
vektion, der Leitung und des Warmeiibergangs zwischen festem und Gaszustand 
erdrtert. Es wird gezeigt, dass es eine stetige Lésung gibt, und eine konvergierende 
Iterationsmethode wird angegeben. 


(Received: March 2, 1959.) 
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Impulsive Rotatory Motion of a Circular Disk 
in a Viscous Fluid’) 


By Ram PrakasH KaNnwaL?), Madison, Wisconsin, USA 


Introduction 


The analysis of rotatory oscillations of an infinitesimally thin circular disk 
in a viscous fluid has been found useful in the theory of oscillation type visco- 
meters. In our attempt [1]%) at finding the solution of the problem both in the 
case of rotatory as well as longitudinal oscillations of the disk we have used 
the spheroidal wave functions of imaginary arguments. The results thus ob- 
tained are not easily amenable to numerical calculations. In the discussions 
of the other workers [2-4], it was assumed that the drag per unit area exerted 
on the disk by the fluid is the same as that which would be exerted on a disk 
of an infinite radius oscillating at the same rate. They have further attempted 
to include the influence of the edge by making corrections of semi-empirical 
type. 

AzpEITIA and NEWELL [5] have extended the above-mentioned theory of 
infinite disk by assuming that the radius of the disk is large compared with 
the boundary layer thickness. Their theory is, though approximate, an im- 
provement over the earlier work inasmuch as the accuracy, as well as the range 


of applicability, of their formulae are well defined and they have checked them - | 


with the available experimental data, 

In this paper we attempt to give an exact analytic solution, amenable to 
numerical work, of the flow when an infinitesimally thin circular disk of 
radius 7) is given an impulsive moment G 6(¢): — d(f) being Drrac’s delta 
function. 


The Equations of Motion and Their Solution 


Employing cylindrical coordinates 7, @, z, and taking into account the 
rotational symmetry of the problem we have u(r, 6, z) = w(r, 6, z) = 0, while 
u(r, 9, 2) satisfies the linearized Navier-Stokes differential equation [6] 


02u 1 Ov 020 v 1 Ov 
Geil > Ott coatia aya meee iF) 


; 1) This research was sponsored by the US. Army under contract no. DA-11-022-ORD-2059 
with the University of Wisconsin. 
*) Mathematics Research Center, University of Wisconsin. 
8) Numbers in brackets refer to References, page 557. 
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where c* = y/o. The symmetry of the problem also suggests that we may 
consider the problem for the region z = 0; the solution for z < 0 can be ob- 
tained by putting —z for z in the following analysis. 
Multiply by e~*’ throughout the equation (1) and integrating from 0 to co 
we obtain 
070 ie OV 07 v oS 
Ort * y Oy * Oe yt or” 2} 
where v is the Laplace transform of v. 
An appropriate particular solution of this equation is 


pee | (E, s) AG é) e-2Vaersle ; (3) 


which can be generalized by integrating it with respect to € from 0 to ov, 
yielding 


~~ 


B= [ AG.) Jr, 8) eo VAPHE ae. 4) 


0 


The stress component y,, is [6] 


The boundary conditions are: 


(1) The couple experienced by the disk is G 6(#), 
(2) po, and therefore #,,=0 for r = 79, where 7 is the radius of the disk. 


Therefore 
—<3H [lan+sés) fie he §) ara 
S 6 
= —A2k [Yah ys Ales) HOD ag = 6, ‘ 
0 
and 
oy Fi jee+s AE, s) fil, 6) dé=0. (7) 

e * 


The solution of these two simultaneous equations for A(é, s) is inde- 
terminate [7a], but we can make it determinate in the following way. The 
solution for the steady case of this problem is known and for that case the 
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relation between the moment G and the shearing stress #9, is [8] 
| SiGe? 


8 aril —v 


? 


= fO0+7 2 % ; | 


(Do heres in 


0 Sia OF? 52237 ae 


We assume that the impulsive couple G d6(¢) produces the stress 


AGES LOL Hie tage | 
Es 8 arr — (9) 
0 A plone eae | 
The conditions [9] are equivalent to 
| a —— , LOR Ft a | 
[Veereae 8) flr 8) dg =} Sant? (10) 
0) ut LOS Tac ae | 


It may be observed that the equation (10,) is the same as the equation (7). 
The solution to these equations is unique and we have [7b] 


A(é, s) 


3Ge_ {fsiné ry iL 
8Saprn EN cose 7 ro} Vee Ee ee 
We have now to verify that this value of A(é,s) is also a solution of the 
equations (6) and (7). As already pointed out equations (7) and (10,) being the 
same we have only to show that the integral equation (6) is satisfied. In order 
to do that we observe that 


(" Me cee r) =z Vir PMG (12) 
Putting the value of A(&,s) in terms of Jg(7). €) in the equation [6] we get 
Jara(%o> E) Jo( Tey <) Zi y2 
dE — il 
i ee $ 3 \a V9 ( : 


Now it is known that [7c] 


Vol. X, 1959 Impulsive Rotatory Motion of a Circular Disk in a Viscous Fluid 555 


Putting A = 1/2, w = 3/2 and » = 2, it is a matter of simple calculation to see 
that the equation (13) is just an identity. Therefore the value of A(é, s) given 
by the equation (11) is an appropriate solution of equations (6) and (7). Sub- 
stituting this value of A(é, s) in the equation (4) we have 


co four 
_ 3Gec aed i g—2Vet s+ sl 
B= VE Tyr 8) 8) ee aE. (15) 
8 | Ter J Ve &2 4 5 
Now 
Fi o+tc0 pst-2Vee rele i / eto pvt-2Vovle 
aay | Py a ae aS 5 Guess 7 —_.— dy, (16) 
G-1t% = cao , Setec Vy 
where 


y=C@l@+s5, 6=CP+.). 


The value of the right-hand side integral of the equation (16) is known and is 
given as [9] 


6 +%OO 
1 ~ yt—zV y/c —22/4c%t 
oS) of = dy => Gi reveely (17) 
LIE if he Vv /x t 


Therefore from (16) we have 


o+10 


4 _xVptee Ic —c® Et 22/4024 
1 est zVc2&+s/¢ é 
fs / aes is a (18) 
cae oe Ge s mt 
rn bok: 


We are now able to invert the equation (15) and obtain 


[o.@) 


eae 27/47 t eas Pe 
3 = agate [Oo VE alto d) Hira 9) 


To evaluate the integral on the right-hand side of the equation (19) we 
__ observe that [10] 


: lo) 2m + 5/2 2 
. 1\5/2 5/9 § : -( ayn 
_ Sr2(to &) Al” §) = (=) " ae m! I (m + 5/2) 


(20) 


3 the: 
x oh (—m, —Z — m2: are 
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Therefore 
2 £2 12 1 \5/2 /2 = Renate goa v5 \™ 
CP ET os Pade) dale £)= (5) % o>, m! I’ (m + 5/2) ( rs an 
s ore 
x ory (—™, cae = Hibs Ze ay 
We are thus required to find the value of the integral 
rs ghee ade, (22) 
0 
which can be readily done by the help of the relation > 
oe ete dé = — 1 f gamer gtr dé. (23) 
0 0 
Thus 
m —tce2&2 m+ 1)! 
fe +3 6 bore aé = ee (24) 
From the relations (24), (21) and (19) we finally have 
SG re ae a m+ 1 ve \m 
ere eas mx pe £512 ¢3 I (m + 5/2) ( thr) 
vd (25) 


x A (—m, — 3 — m2; 4). 


2 
2 Y5 


Since ,F, (— m, — 3/2 — m; 2; r?/r() are polynomials, the value of v can be 
calculated to any desired degree of accuracy. In fact if we put 


. _ 3Yn mt v2 \m 3 7 
In = 4" Pome oy (—aes) A (—m— 3 — meee), 


then we have 


th 75. | i Safe 
= tare a v3 eC aaa ras? 
27) 
hi Zl d Oo mma HOS sears ( 
ee el ie. 
Is a Be Ao pea eT 64 at 


= v8 99 4 231 8 ZS) ae 
= oso tus Fe Se eyes =| 
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Angular Velocity of the Disk 


The angular velocity of the disk is given from the relations (25), (26) and 
(27) as 


gi a 
Yojrar, 32 m2 ut c3 | 20 tc? 1120 22 c4 
2-0 at) 
(2217) 323 v8 \ 
8064 £3 c8 eTse¢ re fe 
This result can be written in a closed form [7d] 
eel G pe as ele r\4 
M0 = 32 SR py Bee 22 (3 te ee | i) 
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Zusammenfassung 


Es wird eine analytische Methode fiir die Bestimmung der Geschwindigkeit der 
Strémung und der Winkelgeschwindigkeit einer kreisformigen Scheibe gegeben, die 
in einer zdhen, inkompressiblen Fliissigkeit einen Drehstoss erhalt. 
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On Relative Measurements of the Viscosity of Gases 
by the Oscillating-Disk Method 


By JosepH Kestin!), WOLFGANG LEIDENFROST 2 andC, Yew) 
Providence, R.I., USA 


1. Introduction 


Until recently it was impossible to perform’absolute determinations of the 
viscosity of fluids by the oscillating-disk method owing to the absence of an 
exact theory which included the edge effect in a rigorous way. Such a theory 
was recently provided by NEWELL [1]*) who obtained a series approximation 
for the case of a thin disk oscillating between two fixed plates. NEWELL’s 
theory can be used for absolute measurements on condition that the boundary 
layer thickness, 6, is not smaller than a certain specified quantity. A thorough 
experimental verification of this theory was undertaken and will be described 
elsewhere [2]; it provided an excellent experimental verification of NEWELL’s 
result. 

It must, however, be realized that the restriction on the magnitude of the 
boundary layer thickness, 6, i.e. actually on the magnitude of the dimension- 
less ratio 6/b, where 6 denotes the separation, makes it unsuitable for an im- 
portant class of conditions which are otherwise desirable from the experimental 
point of view. Consequently, it is still necessary to devote a good deal of 
attention to relative methods of measurement which rely on suitable calibra- 
tion procedures. 

In an earlier paper [3], Kestrn and WANG provided working formulae for 


the relative method, treating two limiting cases, namely those when the 
dimensionless separation 


B ema. (1) 


is either very small or very large (‘infinite’ spacing). In certain applications, 
particularly when operation at elevated temperatures is envisaged, it is found 
necessary to use intermediate separations, and it is useful to extend their 
‘method to include this case as well. Furthermore, with the results of [2] now 
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available it is possible to provide a very precise verification of the dimensional 
argument given in [3]. This showed that the dimensionless ratio 


coon 255 (2) 


of the actual viscous torque M experienced by the disk to the ideal torque 
M,, evaluated without regard to edge effects is a function of the geometrical 
dimensions of the system, all expressed as dimensionless ratios with respect to 
the boundary-layer thickness. Since the values of viscosity in [2] are accurate 
to 0-05%, the test is a very critical one. 

The present paper contains a derivation of the extended formulae and an 
experimental verification of the soundness of the calibration procedure. 


2. Working Formula for Moderate Spacings 


The characteristic equation of the motion of the system was given already 
in [3], but will be repeated here for ease of reference in a slightly modified 
version. It has the form 


A)? + 1+ a a) 3/2 2 [coth (B, st) + coth (B, s1?)] | 
K,(& st!) | ai | (3) 


2 d 
me Ro Ky (& Sue) 


The preceding equation has been written for a disk of thickness d and radius 
R = &,6 oscillating between two plates separated from it by distances , = 8,6 
and b, = f.6 respectively, the lateral separation being ‘infinite’. The factor C 
which can be regarded as a function of 6/b, with 6= (b, + 0,)/2, has been 
defined in equation (2); s denotes the complex frequency of the motion, A) 
is the damping decrement of the wire and K, and K, denote the modified 
Bessel functions of the second kind and of orders 1 and 2 respectively. 
Finally, 6 denotes the boundary layer thickness given by 


s=5|/— , (4) 
ag 
where y is the kinematic viscosity of the fluid, and a, is the natural frequency 
of the oscillating system whose moment of inertia has been denoted by J. 
For large lateral separations [3], we may approximate 
Rue ath ¢ AN ius 


Ky(Eo Sul) g es sil2 ’ 


and the remaining task is to find a suitable approximation for coth(p She), 
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remembering that the complex conjugate frequencies s and s are given by 


where 9 = T/T, is the ratio of the actual period T to its value in a vacuum, 
Tj, and where 


Wee ae (5) 


2am Gn +m 


is the decrement of damping, taken with respect to two amplitudes, «, and 
hn +m, M full cycles apart. 
We note that 


=” [coth (B, st?) + coth (By st/?)] 


sinh (X + Y 1) (6) 
~ cosh (X + Y 2) — cosh[(X + Y 2) (8, — Bs)/(B1 + Bo)] ’ 
where 
1 
X,Y = gow [LF 44+ e4)| +h. (6a) 


It is easy to see that in most realistic installations the difference |6, — B,| is | 


much smaller than the sum f, + fy. Denoting the mean separation by 


B= > (Br + Ba) 


and assuming 


we find that 


1 sinh (X + Y 7) 
Be 1/2 eNews 
[coth (B, st?) + coth (B, st2)] x cash Goda eae 


if we put 


cosh [(X + Yi) eal al (7a) 


We find from equations (6) that 


mee de) s(t baad 


so that 
A+Vi 
Pi + Bs 


= O(1). 


[Po— Bil<<28, (7) | 
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Hence, the error in the approximation (7a) is of the order 


cosh (B, — By) — 1 


compared with the leading term-cosh (X + Y7) in the denominator of equa- 
tion (6). Taking the series expansion 


cosh (8, — Ba) — 1 = 5 (By — Ba)? + py (Br — B+ « 


we prove that the error is quadratic in (8, — B,). Supposing that the disk is 


_ exactly in center at room temperature (say 20°C), and that it moves off center 


so that |B, — B,| becomes as high as 5% of (B, + Bs); which is very unlikely, 
we find that the error will amount to only 0:0025% of the term retained. The 
effect of this on the viscosity will amount to + ie -005°%, because the latter 
is proportional to 6%. In most calculations an error of this magnitude can be 
ignored. 

Substitution into equation (3) and separation into real and imaginary parts 
yields for the imaginary part the following working formula: 


eg — 4) —C 2 |, Ka HK, + “py (H+ 36 zee)|- 0, 


6\6 fs 2R0 
where 
1 — (3/2) 4 — (3/8) 4*@ gy «1 + (3/2) A — (3/8) A t (8) 
A, = 21a Ase ’ 2 21/2 @3/2 ? | 
be sin Y > sinh Y 
Rae cosh X — cosY ’ rae cosh X — cosY * | 


Owing to the complexity of the resulting expressions, calculations must be . 
performed with the aid of a digital computer. 

The formula in equation (8) shows that when moderate separations are 
used, but provided that condition (7) is satisfied, i.e. provided that 


[b.—0,|<26, | (7a) 


j the evaluation of the results demands only the knowledge of the sum of the 
_ separations 25 and not of each of them separately. This is very convenient, 


because it is much easier to measure it accurately. It can be predetermined 
when proper spacers (such as Johanssen blocks) are used, as in [2]. 
3. Further Verification of the Hypothesis 


Using equation (8) and the data in [2], the edge correction factor could be 
_ calculated for the 132 experimental points provided in it. The gases range 
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from the very light hydrogen to the very heavy xenon and include: air, argon, 
carbon dioxide, deuterium, helium, hydrogen, krypton, neon, nitrogen, oxygen, 
water vapor and xenon. The result is shown in Figure 1 as a plot of C against 
the dimensionless boundary layer thickness 6/6. This is a slight departure from 
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Figure 1 
Edge-correction factor C for twelve gases: small separation. 


[3] where the variable x) = 6/R) was used. It is essentially equivalent because 
for a given arrangement the previous variable is proportional to the present 
one and differs from it only by the constant factor b/R,. The diagram in 
Figure 1 fully confirms the hypothesis as not a single one of the 132 points 
deviates from a smooth curve by more than 0-1%. Furthermore, the curve 
passes through the point C= 1, 6/b=0 and is asymptotic to NEWELL’S 
theoretical value 


C aeatOree 


The characteristics of the suspension system used in [2] are repeated here for 
convenience: 


R = 35-021 40-05 mm, 
d = 0-8258 + 0-00038 mm, 
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I = 43-1672 + 0-0063 gcm2, 
b, = 1-0094 -+ 0-00254 mm, 
by = 1-0094 + 0.00254 mm, 
Ty = 25-757 40-001 s. 
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Figure 2 
Edge-correction factor C at two temperatures: moderate separation. 


A further verification is provided in Figure 2 which contains measurements 
performed with nitrogen at 25° C and 75° C at intermediate separations as well 
as one point obtained with helium. The suspension system had the charac- 
teristics: : 
R = 34-994 + 0-0025 mm, 
ad =1.0992 + 0-0005 mm, 
I =56-740 + 0-087 gcm?, 
b, + bg = 7-7980 + 0-005 mm, 
T, = 34-805 + 0-001 s, 


and the values given by MicHELs and Gipson [4] were used in the evaluation. 
- Once again, a single curve is traced at the two temperatures; the curve passes 
through C = 1, 6/b = 0, but is no longer asymptotic to NEWELL’s value 

‘ Cus 1-3059". 


oo 
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This is as expected, because the latter value is applicable to small separations 
only. 

In concluding, the authors wish to note that although operation at large 
separations 0 is desirable, there seems to exist a natural practical limit for it. 
When the separations 0 are increased unduly, in the authors’ experience to 
about 50 mm, the reproducibility deteriorates to an unacceptable degree. This 
seems to be due to convective currents and instabilities caused by small 
temperature differences throughout the mass of the gas and accentuated by 
the evolution of frictional heat in the boundary layer. At large separations the 
Rayleigh number seems to exceed its critical value for stability. 
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Zusammenfassung 


Die Abhandlung enthalt die Ableitung einer Arbeitsgleichung fiir Relativ- 
messung der Zahigkeit, bei der nach der Methode der oszillierenden Scheibe eine 
diinne Scheibe zwischen zwei feststehenden Platten Drehschwingungen ausfithrt. In 
der Formel ist der Randeinfluss nach Kestin und Wane durch den sogenannten 
Randkorrekturfaktor beriicksichtigt. Dieser aus dem Verhidltnis des wirklichen 
Reibungsmomentes zum theoretischen sich ergebende Faktor ist abhangig von der 


Geometrie der experimentellen Anordnung, deren Abmessungen mit der Grenz- 


schichtdicke dimensionslos gemacht sind. | 


Da von Kerstin und Wana fiir Anordnungen mit sehr grossen und sehr engen 
Spalten frither schon Formeln gebracht wurden, beschrankt sich die Ableitung hier | 
auf mittlere Spaltweiten. | 

Um die Richtigkeit der der Arbeitsgleichung zugrunde liegenden hauptsdchlichen 
Annahmen zu uberpriifen, wurden die nach dem Relativverfahren gewonnenen 
W erte mit neuen und sehr genauen kiirzlich von Kestin und LEIDENFROST durch-. 
ae Absolutmessungen verglichen. Es ergab sich eine Ubereinstimmung von 

’ 0° 
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Thermoelastic Problem for an Isotropic Sphere with 
Temperature Dependent Properties’) 


By Jerzy Nowrnskt, Madison, Wis., USA 2) 


Introduction 


Modern structural elements are often subjected to temperature fields of 
such intensity that their mechanical behavior approaches that of viscoelastic 
bodies, and a representation by means of a classical Hookean model does not 
adequately describe their actual properties. 

In spite of this fact, it was shown recently that the elastic thermostress if 
based on temperature dependent properties of a material corresponds to the 
upper value of stress in an inelastic problem for the incompressible viscoelastic 
Maxwell body [1]}%). The results obtained for rotationally symmetrical plane 
states of stress and strain with steady-state radial temperature gradient make 
it plausible that also in other cases the elastic thermostress, defining the initial 
conditions of a viscoelastic thermal process, might serve as an estimate of the 
decisive thermostress in actual structural elements. 

Such a postulate will be also adopted in the sequel for metallic elements 
under consideration. 

On the other hand, experimental data for common steels evidence that at 
elevated temperatures Porsson’s ratio approaches its upper bound for iso- 
tropic materials namely 0-5 [5]. Of course, this is only approximately true and 
the respective values oscillate between 0-4 or 0-45 and 0-5. 

The assumption of the elastic incompressibility of a body at elevated tem- 
peratures, however, permits us to obtain a simple solution to the problem, in 
a closed form convenient for numerical computations. Moreover, the solution 
may pertain to any polarly symmetrical temperature field T = T(r), 7 being 
a position vector, as well as to any variation of Youne’s modulus E = E(T), 
shear modulus G = G(T) and coefficient of thermal expansion « = «(T), with 
temperature. 

This assumption will be, therefore, adopted in subsequent arguments as the 
second basic postulate. 


1) This research was sponsored by the United States Army under Contract No. DA-11-022- 
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As far as the influence of the temperature dependent properties of a 
material on the thermostress is concerned, it has been investigated, to the 
present writer’s knowledge, in the case of an infinite thick-walled tube and a 
circular ring, only in four papers [1-4]. The first two papers treat elastic in- 
compressible material, the last two an arbitrary but constant Poisson ratio*). 

The scope of the present paper is confined to an extension of the previous 
solutions to the case of a solid or hollow sphere subjected to spherically sym- 
metrical temperature fields. Formulae have been derived for thermal stress 
dealing with arbitrary variation of temperature as well as arbitrary variation 
of Younc’s modulus, coefficient of thermal expansion and thermal conduc- 
tivity with temperature. 


Strain and Stress Fields 


Consider an isotropic elastic sphere, with the pole of spherical coordinates 
at the center of the sphere, subjected to spherically symmetrical external 
tractions and spherically symmetrical temperature field T = T(r), with 7 as: 
position vector. Assume that Youne’s modulus E and coefficient of thermal 
expansion of the material of the sphere show marked dependence on tempera- 
ture, and that Porsson’s ratio y is constant and equal to 0-5. 

With usual notation, we obtain, in the case considered, the Hooke- 
Duhamel equations 


whereas all other components of stress and strain tensors vanish. In fact, T(7) 
denotes the excess of temperature above the uniform ambient conditions at 
which thermostress disappears, and elastic and thermal coefficients have their 
usual values. 
With w as radial displacement we have for the only non-vanishing unit 
elongations in the radial and tangent directions, 
ae 2 
¥ ar 


(2) 


Since the elastic dilatation of the body is zero, by hypothesis, equations (1) 


“a Addendum by Proofreading: The author wishes to express his thanks to Professor CHANG 
University of Minnesota, for drawing attention to the paper Stresses in a Metal Tube Under Both 
High Radial Temperature Variation and Internal Pressure by C. C. Cuane and W. H. Cuu, J. appl. 
Mech. 21, 101 (1954), in which a related problem has been investigated on the basis of experimental 
data. Also a recent paper by R. Trostrex, Stationdre Wéarmespannungen mit temperaturabhangigen 
Stoffwerten, Ing.-Archiv 26, 416 (1958), giving a general solution to the thermoelastic problem for 
temperature dependent material properties, using the perturbation method, should be mentioned. 
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and (2) yield simply 
d 
re fat at (3) 


with a notation 
T(r) 
= | a(r) dts (4) 
0 


An integration of the differential equation (3) yields the radial displacement 
3 C 
OS ea (5) 


where C, is a constant of integration to be determined later from the boundary 
conditions, and ‘ 


Nail") = [oo do; (6) 


@ is any convenient lower limit for the integral, such as the inner radius of a 
hollow sphere. 


Clearly now 


oe 3 a 
&, = Sg [3 Xar(”) 35 eal =r C,| , ey — PEE [3 vt) oF Cy] 2 (7) 


and the second equation (1) gives 


i mie Oy) Oa) (8) 


y> 


In addition to the foregoing equation, in the absence of body forces, we have. 
the following well-known stress equation of equilibrium 
do, 2 


dv sae O,) = 0. (9) 


By virtue of these two equations a trivial integration gives finally the radial 
_ and hoop stresses 


o, = 4 [3 Y(r) — wilt) + C1 Lao] +. Ce» (10) 


“a, = 6 [2¥,0) + 2D 7,0] —2 12 v.(0) + BT) _ 
+2[2 neal) +2] C402, 
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with the notation 


Af 2 4 @D 
v= [ DP [erode dc, yl) =/ BD) > de. 
a a (12) 


and C, as constant of integration. 
Thus we have obtained, as was announced earlier, a solution in a closed 


form and valid for any symmetrical temperature distribution as well as for any 


variation of T(r) and a(r). 
Let us now investigate two characteristic cases. 


Solid Sphere 


In this particular case the lower limit a of the integrals may be put equal 
to zero. 

For y= 0, because of symmetry, wo ought to have «= 0. This yields 
C, = 0 in equation (5) since by virtue of the theorem of the mean value of an 


integral, 
lim —y Yo?) = 0, (13) 


with 0 substituted for the subscript @ in y,,. A similar notation will be sub- 
sequently used, for a = 0, in the symbols Y,,, y, and y,5. 

We proceed now to show that o, and o, also remain bounded, for a > 0 
and r +0, if we take C, = 0. Im fact, since, for a = 0, the integrands in W,(7) 
and y,(7) increase infinitely by approaching the lower limit of the interval 
of integration, we have to compute the difference of the improper integrals 
P(r) and y,(r), 

D,, = 3 Fr) — »,(r) (14) 


assuming @ + 0, and then pass to the limit by allowing a >0. This yields’ 
a . 
limD,, = — E(0) (0). (15) 


a—0 4 
10 ; 


wo] 


Now, in order that 7,2 in equation (10) and the last bracketed term in (11)9 
remain bounded as a > 0 and 7 + 0, it suffices to put C, = 0, as was asserted. 

Thus, finally, we obtain for the stress at the center of a solid sphere the 
bounded expression 


(0) = 0,(0) = — 5 E(0) ®0) +c. (16) 


—_—* 
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Let us write explicitly the difference 6, — 6,, for C, = 0 and a= 0, from (8) 
in the form 
[42 QGdg 
05 FO. 2 E(T) Oe > —__——— . (17) 


ey 3 
=~ ar 
=) 


Clearly the second bracketed term in the foregoing equation represents the 
mean unit thermal expansion of the sphere with radius 7. Hence, the follow- 


ing theorem holds: In a solid sphere with temperature dependent properties 
_ subjected to radial temperature gradient, the difference between the local 
_ hoop and radial stresses at a point 7 is in direct proportion to the difference 


between the local value of the thermal expansion and its mean value within 
the sphere with radius r. The factor of proportionality is double the local value 
of YouNG’s modulus. 

Since half a difference between each two principal normal stresses at a 
point designates a principal shear stress at the point, the foregoing theorem 
provides us with direct data for the Tresca-Saint Venant or Huber-Mises- 


_ Hencky criterion of plastic yielding at the point. Furthermore, the factor of 
_ proportionality being equal to the modulus of elasticity, we can at once 


perceive in what degree the temperature dependent properties of a material 
influence its tendency to plastic yielding. It appears, at this early stage of our 
investigation, that the assumption of a temperature-independent elastic modu- 


lus and the use of its value corresponding to the ambient conditions lead to 


an overestimation of the intensity of stress, provoking plastic yielding at 
elevated temperatures. 

The constant of integration C, can be now computed from the boundary 
conditions at the surface of the sphere 7 = b. For a tractionfree surface, for 
instance, we get 


o, = 12 [Po(r) — Yo(5)] — 4 Lol”) — Yol2)] . | 


oy = 12 [%olr) — Py(0)] — 4 [polr) — yol6)) + 2 [2 yo) — £(2) 9]. | 


Sphere with a Central Hole 


Let us denote, in this particular case, the inner and outer radii of the hollow 
sphere by a and 3, respectively. Suppose, furthermore, that in addition to a 


_ steady-state radial temperature gradient, the spherical container is submitted 


to an action of an internal uniform pressure p, and embedded in an elastic 
medium of Winkler type characterized by the modulus %. These assumptions 


clearly yield, o,(a)=—p, 9,(b) =—x u(b) . (19) 
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From (5) and (10) a trivial computation yields the following values of the 
constant of integration: 
3 Fay — Vad ae = as Xa1b os 
C= SOs (20) 
Xa2b e aa 


an additional subscript in the notation of the functions (6) and (12) being 
required to indicate their value on the outer surface 7 = b. 

For x = 0, one deals with an instance of a sphere free from external load on 
the outer surface. 

For x = oo, an instance of a container embedded in an abeohitele rigid 
medium is obtained. In this particular case 


OO eae el (21) 
and u(b) = 0. 

The elastic medium which, in a general case, surrounds the container may 
be, of course, represented by a second hollow sphere, constituting an outer 
layer of a composite double layer sphere. It is easy then to compute the 
respective modulus of the medium x. 

To this purpose, let us assume that the medium is also temperature sensi- 
tive, and denote all quantities concerning the medium by an asterisk. We 
can adapt all the foregoing equations by substituting 6 for a, for the inner 
radius of the outer sphere, and c for b, for the outer radius of that sphere. 
Moreover, we have to take x = 0 and denote inner pressure # by x. This 
procedure leads to the following equations: 


3 6 
uF = a [e 2 @D* do + 


D 


ry 


os 
o* = af ks an 0? B* do toa] Br ) do (22) 


+407 [ =I do+ce 
b 


T(r) 
with @* = ‘i a*(z 


Using Piston similar to that used previously one gets from (20) the 
constants of integration, 


Ch =e = OF ae. (23) 
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Now the condition w*(b) = 1, which expresses the true nature of the modulus 
x, as being a pressure which provokes a unit radial displacement on the inner 
surface r = b of the outer sphere, yields 


n= 4 (3 ae a Voc my Noe b) ‘ (24) 


If the improper integrals in the last equation converge for c -> oo, then the 
limit of the right-hand member in (24) will represent a unit response of a 
temperature sensitive infinite elastic space produced by an expansion of the 
(inner) sphere. It follows that we shall be concerned with a thermo-elastic prob- 
lem for a temperature sensitive container embedded in an ideal temperature 
sensitive elastic medium. 

Let us now take the situation as it often occurs in actuality, namely that 
of a medium at zero excess temperature, 7* = 0. Then ®* = 0, E*(0) = E%, 


and one obtains simply 
4 Ex 


3b” ee 


+= 


which is a known result [6]. 


Temperature Field 


Let us suppose that the thermal conductivity of a material k has a marked 
dependence on local temperature. Then the general equation of the temperature 
field, for a steady-state conduction in a system free of sources and sinks, 


becomes 
RV22 == Vee VE =O; 


with V and 7? denoting a gradient and a Laplace operator, respectively, and 
the dot denoting the scalar product. 
In the polarly symmetrical case we get in polar coordinates simply 


d 2 aT = 27 
— (A ieee ) =), (27) 

or after integration 
frav=-3 +k, (28) 


with K, and K, as constants of integration. _ 
For most practical situations the nonuniformity of k can be expressed 


without any substantial error by a linear function 


k=k,(1—BT), i 
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where k, is the known conductivity at the ambient temperature and f is a 
coefficient, in general, positive for homogeneous metals. 
To be specific, let us assume the isothermal boundary conditions at the 


surfaces of the container discussed, 


y=@, LET, S730, Ae (30) 


Then integration of (28) and solution of the quadratic equation in T yields, by — 
virtue of (30), the following distribution of temperature: 


r-4-V[G-)0-)-G-2)*0-F)- 


Following the procedure of [3] it is easy to show that using TAyLor’s formula | 
for the right-hand member of equation (31) and passing to the limit, with 6 > 0, 
yields the well-known temperature distribution for uniform conductivity: 


r=,+~[n8(1-4)- T,2(1-~)]. (32) 

It will be shown later that for moderate temperature gradients, such as 
assumed in the following numerical example involving prescribed steady 
temperature values on the surfaces of the container, no significant difference 
exists between the temperature fields for temperature dependent and in- 
dependent conductivities®). In view of this, the simple temperature distribution 
(32) instead of the complex one (31) can sometimes be adopted. 


Numerical Example 


Let us consider a hollow sphere under steady-state radial temperature 
gradient and no surface tractions. Using our previous notation we assume that 
a= 50(cm], 6= 100[cm], T, = 500[°C], T, = 400 [°C]. For a steel container 
we take 

“= (1200 + T°C) 10-8 [PC=t], , 
33 
E = (1758.2 — 3-573 (T °C — 400)]_ [kgem-*] . ae 
The last expression has been computed for the interval 750 < T°F < 1000, 


which includes the interval 400 < T°C < 500, from the experimental data 
given in [5] for the S.A.E. 1095 steel. 


°®) It was shown, however, in [3] for a circular tube, that this is not always the case if instead of 


the prescribed boundary temperatures a heat exch Ww dy and the sur 
e ange bet 
unding me di . 8 een a surface of the bo 
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Assume in (29) that ky = 0-11 [cal cm-! st °C} and 6 = 0:0007 [°C-}. 

Then equations (32) and (31) yield the following temperature fields: 
Ti 300 


’ 


as Feateci west sha) (34) 
T, = 1428-57 — |/2,115,012 (1 — "| — 862,242 (1 =f ), 
4 Y 


| < = 
which for specific values of y = 60 and 80 become 


. | rT, | T, Difference 
= per cent 

60 | 466-67 | 465-51 1:16 

80 | 425-00 | 424-07 | 0-93 


In view of an insignificant difference between the respective values of 
temperature just obtained we shall use the simpler equation (32). By substi- 
tuting its right-hand member in the second relation (33) we get finally Younc’s 
modulus in terms of the position vector, 

Ea (2: bier 


v 


) 108. 


Consequently, we compute @ from (4) and then the stresses from (10) and (11). 
They are shown in Figures 1 and 2 being denoted by (temp. dep.). Along with 


0, (temp.dep.) 


Te fgyle (ee e 
o, (¥=z, EE nega: Oe Oacar) 


0,"*(v-4, E-E(0),a= @(0)) 


Figure 1 


The radial stress. 


the curves related to the latter stresses two other pairs of curves are plotted in 
the figures marked by one or two asterisks. Those marked with one asterisk 
correspond to the well-known formulae for a temperature insensitive material, 
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given, for instance, in [7]. There have been computed for Po1sson’s ratio equal 
to 0:5 and for the values of « and E from (33) corresponding to the arithmetic 
mean of the boundary temperatures, that is to 450°C. It may be seen that, at 
least in the particular case considered, the difference between the corresponding 
stress fields is meaningless as far as the accuracy usually required is concerned. 


1000 + Ty” ( v-46 ~Emean £=2 mean) 


ogt*(v-3, E=E(0), a-a(0)) 


0p/ temp.dep.) 


Figure 2 
The hoop stress. 


Hence, in this instance, we could calculate with orthodox equations for tem- 
perature independent material using mean values of « and E. For the sake of 
comparison, in Figures 1 and 2 the stress field for temperature independent 
material has also been shown assuming vy = 1/3 and the values of E and « cor- 
responding to the ambient conditions. Clearly, in this case the actual influence 
of the temperature on the intensity of the stress field would be considerably 
underestimated. Thus, for instance, the maximum compressive hoop stress on 
the inner surface of the sphere would be almost 27% less than that computed 
from the more accurate equation. 

It should be noted that in the foregoing numerical example the temperature 
field is of such an intensity that it induces stress locally exceeding the proportional 
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limit of the material. Consequently, the numerical values of stress in these 
regions may serve only as an illustration of the influence of the temperature 
sensitive properties of the material on the magnitude of thermal stress. 
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Zusammenfassung 


Es werden Warmespannungen infolge polarsymmetrischer stationdrer Tempe- 
raturfelder in einer isotropen Kugel mit temperaturabhangigen Stoffwerten unter- 
sucht. Das Material wird als inkompressibel vorausgesetzt, was einem oberen 
Grenzwert der viskoelastischen Spannungen in einem Maxwellschen Ké6rper ent- 
spricht. Diese Annahme erlaubt, die Lésung in einer geschlossenen Form und fur 
jede Temperaturverteilung und Temperaturabhangigkeit der Materialkonstanten 
und des Warmeausdehnungskoeffizienten aufzustellen. Es werden allgemeine For- 
meln fiir die Spannungs- und Verschiebungsfelder in einer vollen Kugel und in einer 
kugelférmigen Schale gegeben und die Einwirkung der variablen Warmeleitzahl auf 
das Temperaturfeld untersucht. Ein numerisches Beispiel fiir lineare Temperatur- 
abhangigkeit des Elastizitatsmoduls und des Warmeausdehnungskoeffizienten wird 


berechnet. 


* (Received: April 22, 1959.) 
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Eine Kuppel gleicher Festigkeit 


Von WERNER ISSLER, Ziirich 


Einleitung 


Um den Membranspannungszustand einer beliebigen Schale zu beschreiben, 
geniigt es, fiir jeden Punkt der Schale die Richtungen und die Betrage der Haupt- 
spannungen o, und o, anzugeben. 

Spezieli spricht man von einer Schale gleicher Festigkeit, wenn fiir jeden Punkt 
der Schale die Gleichung 

Oy = 09 =. 96 (1) 


gilt, wobei mit o, die Fliessgrenze (bzw. Bruchgrenze) bei reinem Zug oder Druck 
bezeichnet ist. Die Festigkeit des Materials einer solchen Schale ist offensichtlich 
optimal ausgeniitzt, wenn fiir das Fliessen (bzw. fiir den Bruch) die grésste auftre- 
tende Normalspannung massgebend ist. Geniigt hingegen das Material einer andern 
Fliess- bzw. Bruchhypothese, so ist die obige Definition der gleichen Festigkeit 
weniger sinnvoll. 

ZIEGLER?) hat daher kiirzlich den Begriff der gleichen Festigkeit in einer neuen 
Form verwendet. Er ersetzt die Forderung (1) durch die Fliessbedingung von 
Tresca (bzw. Bruchhypothese konstanter Schubspannung) 


Max (04, |g, |¢; — o2|) = % (2) 


oder durch die Fliessbedingung von v. Mises (bzw. Bruchhypothese konstanter 
Gestaltanderungsarbeit) 


ot os — 0; 6, = Gr . (3) 


In der Spannungsebene (0, 0.) (Figur 1) wird die Bedingung (1) durch den 
Punkt A, (2) durch das Sechseck von TrEsca und (3) durch die umgeschriebene 
Ellipse dargestellt. 


Die spharische Kuppel gleicher Festigkeit 


4 Wir betrachten die spharische Kuppel, die nur durch ihr Eigengewicht belastet 
ist. 
_ Wie ZIEGLER?) gezeigt hat, existiert eine solche Kuppel gleicher Festigkeit im 
Sinne der Bedingung von Tresca. Dieselbe Kuppel lasst sich aber auch als Kuppel 
gleicher Festigkeit im Sinne der Bedingung von v. Misxs ausbilden. Die Lésung 
dieser Aufgabe ist im folgenden kurz dargestellt. 
Als Koordinatennetz auf der Kugel wahlen wir das geographische, wobei wir mit 
# die Lange und mit » das Komplement der Breite bezeichnen. Fiir die gewahlte 
rotationssymmetrische Belastung hat die eine Hauptspannung, zum Beispiel oe 
iiberall die Richtung der Meridiantangente und damit o, die Richtung der Breiten- 
kreistangente. Ist a der Radius, ¢(p) die Wandstarke und y das spezifische Gewicht 


1) H. Zrecwer, Kuppeln gleicher Festigkeit, Ing.-Archiv 26, 378 (1958). 
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der Schale, so erhalt man durch Spezialisierung der bekannten Membrangleichun- 
gen) auf unser Problem die Beziehungen 


d : 
do % Sing) — ta, cosp+aytsinyp=0, o,+0,+aycosp=0. 


Zusammen mit der Bedingung (3) 
of + 03 — 6, 0,= 03 


und der Schalenstarke f) im Scheitel bestimmen sie den Verlauf der Wandstarke 
#(p). Man erhalt zwei Losungen, da man fiir den Spannungspunkt (0, o,), ausgehend 
von A (Figur 1), den Ellipsenast AB oder AC vorschreiben kann. 


Figur 1 


Fliessfigur (oder Bruchfigur). 


Lésung 1 


_ Der Spannungspunkt (¢,, o,) wandert auf der Ellipse (Figur 1) von 4 nach B. 
Die Lésung 


cos 2/6 3. 
Beh bq E= (p+ a 


bleibt auf das Intervall 0 < y < 2/3 beschrankt, da ¢/t) fir p> 2/3 gegen co 
strebt. 
Losung 2 


Der Spannungspunkt (o,, ¢,) wandert auf der Ellipse (Figur i) von A nach C, 
Die:Lésung 


cosz/6 73 
Hg) = fo Ee (p — a4 


gilt fiir das Intervall 0 < py < 22/3. 


2) W. Fiucce, Statik und Dynamik der Schalen, 2. Aufl. (Springer, Berlin 1957), S. 26. 
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Interessant ist der Vergleich dieser Losungen 1 und 2 (Figur 2) mit den entspre- 
chenden 1’ und 2’ (Figur 2), die Z1EGLER?) fiir die spharische Kuppel gleicher Festig- 
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Figur 2 
Verlauf der Wandstarke. 


keit im Sinne der Bedingung von Tresca ermittelt hat. Die Lésung 2 ist offen- 
sichtlich Gkonomischer als die tibrigen drei, da sie fiir alle zulassigen Winke .g mit 
der geringsten Wandstarke auskommt. 


Summary 


Thin shells of optimum strength are usually defined by the condition that the 
two principal stresses in any point are equal to the yield stress. This definition is too 
restrictive; it is sufficient that principal stresses satisfy the yield condition of 
v. Mises or the one of TREScA. 

In this paper a proof is supplied for the existence of spherical cupolas of equal 
strength on the basis of v. Misxs’ yield condition, and an exact solution is presented 
for the distribution of thickness. 


(Eingegangen: 12. Mai 1959.) 
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Propagation of Elastic Waves in an Infinite Plate 
of Cylindrically Aelotropic Material 


By Maninpra Mirra, Calcutta, India 1) 


Introduction 


The oscillations of an infinite, isotropic elastic plate have been studied by 
RaYLEIGH [1]?), Lams [2], Saré [3] and others. In this paper free oscillations of an 
infinite plate of cylindrical aelotropy are considered for motions symmetrical about 
an axis perpendicular to the plate. The frequency equations obtained were used to 
calculate numerically, for the first two modes of propagation only, the dependence 
of the wave-velocity on the wave-length of progressive waves in such plates. 


Solution of the Problem 


Cylindrical co-ordinates are chosen in such a way that the free surfaces of the 
plates have equations z = + a in the equilibrium position. It is assumed that the 
oscillations are symmetrical about the z-axis and that the stress-strain relations are 
given by 


WW = Cyy Cy + Cyp C99 4- C13 Eze » 

68 = C45 Crp + C13 C00 + Cra Exe » | x 
ZZ = C43 Cyy + Cya C69 + Cas Cze » | w 
OZ = Cyy lpn, = Cay ly, VO= <u *2 eg, 


where rv, 00, zz etc. are the components of the stress and ¢@,,., ég9, €,, etc. the compo- 
nents of the strain. Using the well-known expressions for the components of the 
strain in terms of the components (uw, v, w) of the displacement, noting that for » 
motions symmetrical about the z-axis v = 0, and substituting in the equations of 
motion, one obtains 

- 0?u 0?w 0?u 
[en Set tS + cua ee | + Cae + ead Benet oe aie | 


(2) 
O2u 1 | vs O2w 0 “| 02w 


Caa w gece 
(1s + Cas) (5 dz | ry Oz (44 Oy? Site ype ie Bea COR 


the other equation of motion being satisfied identically. 
On putting 
Mia Ae TAG r) effet tevez w = B’ Tlé r) erSctteye 


1) Department of Mathematics, Jadavpur University. 
2) Numbers in brackets refer to References, page 583. 
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A’, B’ being constants, it is easily seen that the differential equations (2) are satis- 
fied if ; 
(Cas V2 + @ C2 — Cy) A’ — y (Crs + Cag) B= 9, (3) 
(C43 + Cas) A’ + (C53 y2 + 0 0? — Cyy) B= 0, 
that is, if 
(Cag 2 + @ 6? — C41) (Cog v? + @ 6? — Cas) + Y? (Crs + Cas)? = 0- (4) 


Let +h,, +k, be the roots of this equation in y. The general solution for (#, w) 
may then be taken as 


pee [A ef haz Je B ea sz ot C eb? a D oe Kié r) eset ; 


5 [“ i+ ee re Pee {A oot _ B gts 
hy (C13 + a4) (5) 


C4a ki Se os — Cy, {C eS 2 SP eta To(é r) eset , 
hy (C13 + Cqa) 


the expression for w being determined with the help of equation (3). 

The boundary conditions at z = + a are the vanishing of the stress components 
22, 2v, 20. These conditions, written down with the help of (1), (5) and the expressions 
for @,,., €gg, etc. in terms of (u, v, w), are satisfied if 


Cag hi + oc? — Cy : 
[i hy (Cys + Caa) emma 


(6) 
Cy hk? +002 —C : 
pi fee ce 1. C, sinh haeeo? 
| ‘ Ry (C13 + Cas) lem) 
jp hs 
(7) 
CrRs 2 
| [ers Ota ee al C, cosh(é ky a) = 0, 
13 44 
_ Sali + Qe — by 
pcan eG 
(8) 
Cy hi + oc? —c 
a k 44 V1 =| 1D h ii 
[a ay tt] Di coshlé A, a) = 0, 
yr, Cag hf + oc? — : 
[a +4088 5 a = qq al = ak hs y) 
9) 
Ople pee ( 
i ee Poe ree “| D, sinh(ék, a) = 0, 
13 1 4g 
where i 
ee Ea St ae Cue Es Coe 
2 ? B ) eee. tba? ee aa D = ip mee 


The above equations will be satisfied fo diff Re 
following conditions hold. T cilterent values of 9 c? if one of the 
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(i)-B, =D, = 0 abd 


tanh(aék;) ky ashi -eck+c, 
tanh(a & h,) hy Cyg RF — OC? + Cy | 
10 
Cys (Cys ot tag) = C41 Cg3 + Cg Sag RZ + Cog QC? ee 
C43 (C13 + C4a) — Cyy Cg + Cygg Cag AZ + C33 9 C* * 
In this case 
u cc [A, cosh(§ h, 2) + C, cosh(é Rj 2)] 
and ; 
Caaf? + 0 6 2— EO oa Cy hk? + oc? —c 4 7] 
Wo . tA, sinh(é h, 2) + —# — au CG, sinh(é Ay 2)| ; 
Hence U,_.=U,__, and nt pang Aa 
Thus the oscillations are symmetrical about the plane z = 0. 
(ii) A, = C, = 0 and 
_tanh(a § hy) ky eis hi— eC + cy , 
tanh(a & k,) hy. 643 8? — 9.6? +6.) | 
(11) 
C43 (C13 + Cas) — C11 C33 + C33 Cag RZ + Cag QC? | 
C13 (C13 + C4a) — Cy C33 + C33 Cag KZ + Cy3 QC?” 


It is easily seen that this frequency equation corresponds to vibrations anti- 


symmetrical about the plane z= 0. 
If h?, say, is negative, on putting h} = — h;?, and writing 


tanh(#,@&é)=ctan(kiaé), hy=whi, 


one obtains equations (10) and (11) in terms of real quantities. 
Similar changes are made if k? is negative. 
As a —> ov, the above equations reduce to 


2 2 2 2 
ky Las Wji-—eP+ty Cs (Cis + C44) — C11 C33 + Cas Can RE + Cys QC 
2 2 2 
hy Cyg hf — OC? + Cyy C43 (C13 + Cga) — 11 C33 + C5 Cag ht + Cy3 QC? 


which is the equation giving the Rayleigh wave velocity in the medium. It is as- 


sumed here that h? > 0, kh? > 0. 
These equations reduce to the corresponding equations in the case of isotropy if 


Cy = beg =A 2M, Cy =lyg=A, Cu=p- 


The oscillations considered correspond to periods 2 2/&c, and 2 2/& is approxi- 
mately the length of the wave-form of the surface of the plate for large values of 7 
since the asymptotic expansion of J,(&7) and J,(€7) involve cos( 7) and sin(& 7). 
However, if the propagation of elastic waves are considered, 2 7/§ is the wave-length 
and c the wave-velocity for large values of vy. The dependence of ¢ on & has been 
determined numerically from equations (10) and (11) for the first two modes of 
propagation only. It has been assumed that 


C4y = 2694, Cyg = 2-363, yg = 9-961, Cg = 0-661, Og = 0°653. 
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Figure 2 
Anti-symmetrical vibrations. 


Vol. X, 1959 Kurze Mitteilungen — Brief Reports — Communications braves 583 
These are such that the values of h? and k? obtained from equation (4) are 
always real. (For transversely isotropic materials these values correspond to beryl.) 


The values of Ve ¢ are given with corresponding values of a & in the following 
table, and plotted in Figures 1 and 2. 


Symmetrical case Anti-symmetrical case 
ag Voc a€& Vo c 
fone) 0-7724 0 0 
5-369 0-8081 0-636 0:5477 
2-616 1-0000 1-394 0:6325 
1-855 1-2247 2-407 0-7071 
1-463 1-4142 6-457 0-7681 i 
1-246 1-5000 foe) 0:7724 
0-990 1:5492 13-263 0-8367 | 
0-361 1-5811 5-806 1-0000 
0 1-584 3-802 '1-2247 
3-260 1-5906 2-926 1-4142 
2°951 1-6125 2:254 1:5811 
1-515 1-7321 1-647 sey CVA 
1-289 1-8708 1-277 1:8708 
1-174 2-0000 1-064 2-0000 
1-005 3°1623 0-473 3:1623 
0-644 5:-4772 0-307 5-4772 


A comparison with the corresponding curves in the isotropic case obtained by 
previous investigators shows that the curves are similar in the two cases (see 
reference [4]). Thus, the propagation of elastic waves in isotropic and cylindrically 
aelotropic plates are similar in character. 

The author is grateful to Prof. B. B. SEN for his kind help and guidance. 
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Zusammenfassung 


Die Abhangigkeit der Wellengeschwindigkeit von der Wellenlange elastischer 
Wellen, die sich in einer unendlichen Platte mit zylindrischer Aelotropie ausbreiten, 
ist auf numerischem Weg aus der vorgangig hergeleiteten « Frequenzgleichung » 


ermittelt worden. Die Ergebnisse sind denen des isotropen Falles ahnlich. 
: Ene 
(Received: May 22, 1959.) ‘ 
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Constant Shear Flow Past Two Circular Cylinders 


By Laxsumt N. Nicam, Kharagpur, India?) 


1. Introduction 


The problem of two-dimensional flow with constant shear past cylinders having 
various cross sections has received considerable attention in recent times. The 
main interest in the problem lies in studying the location of the stagnation points 
and their shift when a non-dimensional parameter N = U/wc, where U and ¢ are 
a typical velocity and length and is the vorticity, is varied from 0 to oo. The © 
pro blem taken up in this paper has been studied by Ray?) who used the method of 
images to obtain an approximate solution. His solution does not appear to be © 
suitable, however, for studying the location of the stagnation points. In this paper 
the same problem has been studied using orthogonal curvilinear di-polar co-ordinates 
and the location of the stagnation points for different values of N has been found. 
The forces on the cylinders are calculated. 


2. Equations of Motion 


In the absence of external forces, the equations of two-dimensional steady — 
motion of an inviscid incompressible fluid are 


Ou Ou 1 op 
mo dae oy o 0%? fe 
Ue ee PKS 
edhe, ss male Se 0. oy’ G2) 
The equation of continuity is 
Ou ov 0 
ron? Rapa a 


where u,v are the velocity components, p is the pressure and @ the density. In 
terms of the stream function y the equations (2:1) — (2.3) reduce to 


Oy, Vey) _ 
CAE a hs OF (2.4) 
ae u = — dy/dy, v = Ow/dx and * = 02/dx? + 02/dy2. The general solution of 
.4) is 
Vey = f(y) . (2.5) 


In the case when vorticity is constant everywhere (2.5) reduces to 


int shag (2.6) 


1) Indian Institute of Technology. 
*) M. Ray, Bull. Cal. math. Soc. 29, 125 (1937), 
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The undisturbed flow is given by 


a 


Vo LS a Hi (2.7) 


Putting y = yo + yy, the equation for p, is 72 y, = 0. The equations (2.1) and (2.2) 
can be integrated to find # as 


P 


@ 


tro} 


(w+ v2)—woyth, (2.8) 
where & is a constant. 


3. Two Cylinders with Line of Centres Perpendicular to 
Undisturbed Stream 
Two cylinders are placed in the stream (2.7) as shown in Figure 1. The cylinders 
are defined by 
fe = oe, (3.1) 
where z=ctan¢/2, z=*+iv, €=&+i7. The solution of (2.6) which satisfies 
the boundary conditions 


ow 
tt = 0, Bue 
OE Ina (3.2) 
Oy 
bee = 0; 3.3 
0& |n=—B ore 
where 
ee cos& + coshy 
c 
is 
1 ea 
p= Uy— ied y? + oS cosn & (A, coshw 7 + B, sinhn 7) . (3.4) 
n=1 


Making use of the boundary conditions (3.2) and (3.3), (3.4) gives 


< . . . [oe e} 

U csinhe sin€é  @ esinh*s sing. DS. sine etecaaa ae Ease at 
(cosé + cosh«)? (cosé + cosh)? re 
and ; 

U csinhf sin€& @ c® sinh? f sin& 


co 
>) — nsinn & (A, coshn B — B, sinha 8) . 
1 


(cos§ + coshp)? * (cosé + coshp)® — (3.6) 


The constants A, and B, may be evaluated by multiplying (3.5) and (3.6) by 
sinn € and integrating both sides with respect to & between a a Se 
Making use of the integrals (for proof see appendix) 


* — cosn & dé el) a (n sinha + cosh «) 
(cos& + cosh)? sinh? « 
and 
cosm Ed > (—1)" 2 e-"¢ [(m® —1) sinh?« + 3 cosh?a + 3” sinh « cosh «] 


(cosé+ cosha)® = 2 sinh 
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in integration, the equations for A, and B, are 


w cr e-r& m (—1)” 


(m sinha + cosha) 


4 \n-Ek —na 
att ale ee sinh « (3.7) 
=n (A, coshn « + B, sinhn «) | 
and 
In (—1)"*) U ce"? — 2 ee a (n sinh B + cosh f) 
Z ) i = sinh B (3.8) 
= —n(A, coshn Bp — B, sinhn £) . 
Substituting for A, and B,, from (3.7) and (3.8) in (3.4), y becomes 
_  Uesinhy il @ c? sinh? 7 
 “cosé + coshn 2° (cosé + coshy)? 
o 2 (—1)"* Ue | n(y—a) =n (n +B) 
——— _ h 
ae ye cosn & |i Car) (e coshn Bp — e cosh «) 
See (3.9) 
(—1)” oe 


lea me) fn e" —% sinhn B + ne") sinhn « 


+ e~"* cotha sinhn (y + B) — e~”* cothf sinhn (n — x)}]. 


It can be shown that the velocity components due to disturbance vanish at 
infinity (7 = 0, § = m). The forces per unit length of the cylinders are given by 


X=—G peosd as, Y= psindas, (3.10) 


where @ is the angle with the normal makes with the x-axis and ds is an element of 
the contour of the cylinder. On the surface of the cylinders, 


(0? + U8) = (1 sey : (3.11) 


Substituting for p, from (2.8) and (3.11), in (3.10) and integrating, the forces on 
» = « are obtained to be | 


Sot Vi a Pee aap + 3" 2)— 78° cosh * 
; 2 Ft eh 7 cosha 2&5 ait dD" ap Ania), (3.12) 
ng 


@C 
Ce 
2sinh?q ’ 
_ (=1)" we : 
ans area (a +B) [2 n* sinhn B + n e”® cotha — ne~”8 cothf . 


en e-"* sinhn (« + B) Loan —coenn P). 


sinh? 
The forces on 7 = —f are ne OG 


[oe) 
xX#0, ie ea m QC ~ 
> cote + —j- coshB (2050, + 37d, Bees), same 
af 


Vol. X, 1959 Kurze Mitteilungen — Brief Reports — Communications bréves 587 


where 
ae 2 
eo 2sinh?p? 
(—1)" we 
b, = = 
sinhn {a + p 


[2 n® sinhn a +  e"* cothp — ne~"4 cotha 


+ - 


e~"? sinhn (a + B) 4n U coshn « 
sinh? 6 33 wc J 


Stagnation points on the cylinder 7 = a. 

The normal velocity is zero on the cylinder. The stagnation points can be 
obtained by equating to zero the tangential component of velocity on the cylinder, 
ie. 

Ow 
dn 


h 0 


n=a 


This gives 


F(é = N (1+ cos§cosha) sinha (1 + cosé cosha) 


(cos € + cosh)? (cos§ + cosha)§ 


ne (3.14) 
+ D(C, N + D,) COS S10 
n=1 
where N = U/m c is a non-dimensional parameter, and 


mae ert 3) ~n(a+6) 

C,, alia ke a {coshn B + e coshn a} , 
n (—1)” 

Pr sinh (« + £) 


{nsinhn B — ne~"@**) sinhn «a + e~”* cotha coshx (« + f) 


— e*? cothp}. 


The zeros of F(£) for different values of N can be found graphically. The calculation 
for « = log2, f = log3 is shown in Table 1. 


Table 1 


100° 120° 140° 160° 180° 


80° 


60° 


40° 


| E 0° | 20° 


N =0 |—0-18 |—0-47 |—0-56 |—0-73 |—0-93 |—1-55 |—2-24 |—2-08 | 6:58 | 21-11 
N-=1 0-97 | 0-68 | 0-61 | 0-48 | 0-23 |—0-22 |—1-04 |-1:94 | 2:50} 12-81 
N=5 5-58 | 5°30 | 5-32 | 5-34 | 5-33 | 75-09 | 3-77 |—1-38 |—13-16|— 20-45 
eeoolseti5 ie ieloe le bt 7 | 1-22) |e 1-27 e334) le ZOF m0 L459 lg 8.32 


Value of F(€) 


Positions of the stagnation points on 7 = « are given. by 
(1) = 4 149° forN=0, 

)é= + 91°, 4 153° forN=1, 

Ve ech 137° for NS. 

4) &€= + 143° for N = oo. 


(2 
(3 
( 
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In Figure 1 four stagnation points Sj, SiS g ANCES, 


if 


S$ 
eee et O}e: Mb 
‘CH 
D' 


“Figure 1 


i 


Shear flow perpendicular to line of centres. ~ 


are shown for N = 1. As N tends to oo, S,, S, shift towards A, B, ends of the 
diameter A B, respectively. S;, S, move away from the diameter A B until they 
coincide with C on the y-axis. for.certain positive value.of N >.1, whence the 
combined point moves on the y-axis to infinity. 


ae ee 


4. Two Cylinders with Line of Centres Dicer to Undisturbed Stream 
‘q Ge f ue ec 


Two cylinders (Figure 2) defined by (3. 1) are place n the stream) «Ani: “ 
1 

Yo =Vx—- a wx. (4.1) 

we: nes sats Crip od 2 ‘ az Ya 


The stream function for this case can be determined as explained in part 3. It is 
found to be 


Vecsiné feyhent. wm c* sin? & 
(cos € + cosh) 2° (cosé + coshy)? 


y= 


\@ Zs csinn & a ‘ : 
= x? ee (a +p) fe”) sinhn B + e-" +8) sinhn a] 


[e"-) sinhn B + ne-"+®) sinhn « 


ioe 1)"*1 wm c2 cosn — 
+ 
: ee (« + B) 


+ e-"® coth sinhn (n — a) —e~ "* cotha sinhn (n + B)]. 


The forces on the cylinder 7 = « are 


pp Recs ee 
= 7 — Sinha | 2 ag by + SY” (dq Ones — by nya) |, (4.3) 
1 
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: vieh és, 
a 2 Lats 
Vase 2a% + ae ai|- 73 a Feo coshe f2avas Aye, + x Aq Ania + b bea): (4.4) 
where 
mM Cc 
4 = >> 
2 sinh? « 


@ c (—1)"t1 


eC a e"@ sinh n (a + 
xn ~ “sinhn (a + B) re Sem 


sinh? 
+ ne"B coth F 


[2 n® sinhn B — ne”? cotha — - 


_. 4n(-1)" V sinhz Bg 
Mee, Sinkm (a+ B) 


The forces on 7 = —f are 


wt O 


X= —=)* sinhp vt aj bj + ~ (4p Onis — Oh Gnas) |, (4.5) 


oO 
~y 


° cosh B 2 ay ay+ (Gy, Gia + bf ra) , (4.6) 


c’ co 
Y= sal ak ay? + poe ain 
1 1 
where 
hee ae 
.  26mhep” 


, 


* w ¢ (—1)" 
“n ~ “sinhn (a + B) 


tee oe 


: sf : ee \ be . ’ —nB : i i 
xX [2 wn? sinhn « — e”% ates 4+ ne-™4 cotha — — sinh» (« + B) 
fee ; ws sinh? 6 ? 


Beh 48; * iz 4n(—1)"4 V sinhn « 
4 sinhu («+ fp) * 


The stagnation points on the cylinder 4 = —fB may be found by equating to 
_ zero the tangential velocity on the surface of the cylinder. This gives 


F(é) N siné sinhB : sin? € sinh B é 

(cos + cosh)? (cos+ coshB)® | 

a (4.7) 
+ DT (N Ay sinn & + B,, cosn )=0, | 
T 
where : 
= ay _ ae OH dp sinhy B — sinhn a) , 
wo ; 
_4)nH1 

pe (cts 


n~ ‘sinhn ( + 8) 
x [n e~" (4+) sinhn B — n sinha aw + e~"8 cothB coshn (« + B) — e~”* cotha] . 
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The zeros of F(&) for some values of N have been obtained graphically and the 
values of F(é) for « = log 2, B = log3 are shown in Table 2. 


Table 2 


Value of F(&) Value of F(&) 
é |v =0|w=12| v=1| N=5 |N=co 


g |N=0 | v=1/2[N =1]N=5|N=00 


0°} 0-026] 0-026 | 0-026 | 0-026 | 0-000 0°| 0-026] 0-026] 0-026} 0-026} 0-000 
20°} 0-018] 0-048] 0-078] 0-318] 0-060}— 20°] 0-018 |—0-011 |- 0-041 |— 0-281 |- 0-060 
40°|—0-009} 0-061} 0-132 | 0-700 | 0-142 |— 40°|—0-009 |- 0-081 |-0-152|- 0-719} 9-142 
60° |— 0-078} 0-059] 0-197 | 1-296 | 0-274 |— 60°} 0-078 |- 0-215 |- 0-353 |- 1-452 |- 0-275 
80°|—0-224} 0-025 | 0-275 | 2:273|0-499|— 80°|- 0-224 | 0-474 |-0-742|-— 2-723 |- 0-499 
100° |— 0-497 |— 0-064 | 0-368 | 3-831 | 0-865 |— 100° |- 0-497 |- 0-929 |- 1-363 |- 4-826 |- 0-866 
120° |— 0-873 |— 0-181 | 0-512 | 6-055 | 1-385 |— 120 ° |- 0-874 |- 1-566 |-2-259 |-_ 7-803 |- 1-386 
140° |— 1-006 |— 0-077 | 0-851 | 8-278 | 1-856 |— 140 ° |- 1-006 |— 1-935 |— 2-863 |- 10-290 » 1-856 
160°|—0-202} 0-587] 1-376 | 7-693 | 1-579 |— 160 ° |- 0-202 |- 0-992 |- 1-781 |— 8-097 |- 1-579 
180°} 0-657] 0-657 | 0-657 | 0-657 | 0-000 |—180°| 0-657] 0-657] 0-657] 0-657] 0-000 


The stagnation points on the cylinder 7 = —f are given by 


(1) €= + 36°, + 163° for N=0, 
(2) €= 88°, 144°, — 15°, —.1471° -for N= #2); 
(3) = — 9°, — 175° for N = 1, 
(4) § = — 2°, — P78" forN = 5, 
(5) € = 0°, — 180° for N = oo. 


For N= 0 the stagnation points Sj, Sj, S; and S; are shown in Figure 2. 
As N-> oo, points S;, S; tend to C’, D’ respectively. The points S3, S; move 
away from the y-axis until they coincide with B’ for certain positive value of N. 
a greater values of N the combined point moves along the line A’ B’ produced 

0 co, 


Figure 2 
Shear flow parallel to line of centres, 
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In conclusion I should like to thank Dr. S.D. Nicam for his help in the preparation 


_ of this paper. 


s 


APPENDIX 


4 9 


cosn 6 
gif — 
. / cos 8 + cosh ae. 
0 


The recurrence formula is 
Faas + areosho T+ Buy = 0. (1) 


J, and J, can be evaluated easily and their values are 


a6 dd cos 6 dé 
mee : ag ty i; 
, / cos@+coshp — sinhg’ Ny vf (cos6 + coshg) (1 — cothg) . 
a 6 


in (1) put 7, = x", so 

x*4+ 2coshog¥+1=0, 
which gives 

«% = —coshg + sinhg. 


Most general solution of J, is 


I, = (—1)"[4 coshn » + Bsinhz g¢], (2) 
with initial conditions 
ey 
= 3 
In sinh ) ( ) 
I, =a (1 —cothg). (4) 


Constants 4 and B are evaluated from (2), (3) and (4). Substitution gives 
mie 


n 
dymndt#) sinhg ° 


Differentiating it with respect to ¢, one gets 


- — cosn 6 dé (-1)" xe"? 
a4 (cos 6 + cosh@)? sinh’ p 


(n sinhg + cosh@) . 


Again differentiation with respect to p gives 


ns n SP 
cosnOd0 (—1)" ae 
a (cos 8 + cosh q)? 2 sinh® p 


x [(m? — 1) sinh?p + 3 cosh?g + 3” sinhg cosh@q] . 
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Zusammenfassung 


Es wird die Scherstrémung um zwei Kreiszylinder behandelt unter Verwendung 
orthogonaler krummliniger Bipolar-Koordinaten, wobei die Hauptstromung sowohl 
normal als auch parallel zu den Zylinderachsen verlauft. Die Krafte auf die Zylinder 
werden berechnet, ebenso die Lage der Staupunkte in Abhangigkeit von N = U/c. 


(Received: January 12, 1959.) 


Supersonic Viscous Flow Past a Convex Corner 


By Ajit Kumar Roy, Bangalore, India?) 


Introduction 


In the case of supersonic viscous flow past a convex corner the local behaviour 
of the emitted wave near the corner is a pure rarefaction, in which the pressure falls 
monotonically (and for which, aft of the corner, the flow cannot separate from the 
boundary). This is borne out by all Schlieren photographs of PRANDTL-MEYER’S [1]?) 
expansion near the convex corners. Since the corner point lies embedded in the 
boundary layer, a simple mathematical model has been introduced on the assump- 
tion that the effect of viscosity within the boundary layer is as if to round off the 
corner. Figure la gives the picture of the usual Prandtl-Meyer expansion at a 
convex corner in two-dimensional, steady, non-viscous, supersonic flow. Figure 16 
gives the corresponding picture in the viscous case. Figure 1c is the mathematical 
model posed for the present problem and obtained by rounding the corner by 
means of the Prandtl-Meyer stream line which approximates to an arc of a circle 
in the case of a weak expansion (i. e. for very small values of the angle g). Under 
these circumstances, the flow is a boundary-layer flow along a curved boundary 


Figure 1 


a Department of Aeronautical Engineering, Indian Institute of Science. 
2) Numbers in brackets refer to References, page 602. 
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with favourable pressure gradient (and therefore without separation). The usual 
Polhausen method with different assumed velocity distributions is employed to 
determine the flow. 


c* 


Cc 


Notations 


critical velocity of sound; 
constant in the temperature viscosity relation 


—_|/ tw , To +S (assumed 1-5 for the numerical evaluation in the 
Dee Ty +S present case under zero heat-transfer condition) ; 

specific heat at constant pressure; 

enthalpy; 

thermal conductivity ; 

Mach number; 

static pressure; 

PRANDTL’S number; 

absolute viscosity ; 

density ; 


BQ; 


(vy, 6) polar coordinates of a field point; 


Re 


S 


~— 


ROMER TE SSH 


oe og 


part 


gas constant; 

free stream Reynolds number (2 »,9/c*) ; 

SOUTHERLAND’S constant; 

absolute temperature; 

velocity components in the radial and cross radial directions respectively ; 
angle contained between the terminating Mach lines; 

angle contained between the directions of the two walls; 

sheering stress; 

ratio of the specific heats of the air (i. €., c,/Cc,) : 

second viscosity coefficient (given by MAXWELv’s relation 34+ 2 = 0); 
length O, O, (Figure 1d) ; 

length O, O, (Figure 1d); 

length O, O, (Figure 1d) = (i — Dd). 


Subscripts 


free stream conditions (undisturbed); 
local free stream conditions (i. e., at the outer edge of the boundary layer) ; 


conditions on the boundary wall; 
standard conditions of the flow parameters attainable on the first Mach line 


(sonic conditions for Ma, = 1); 
conditions upstream of the flow deflection ; 
conditions downstream of the flow deflection. 


Analytic Development 


The usual stress and strain relations [2] in polar coordinates (Figure 1c) are 
given in the following: 


Strain components: 


Ou. pels ye U . laren ov ay (1) 
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Divergence: 


Ou Ou 
4 Sars 2 
| are Sr oF a GMeae (2) 
Stress components: ‘ ; 
Ou u 1 uU uU 
bp = —PtRAt He, =P tala +7 +7 og) HOH) 


Ou uU 1 Ou 2 ov U 
pope Spa en p+ alse + 4 Se) + a(S spt 25), (4) 


1 ou ou uv 
Pro = Here =H (5 00 F Or aaa 


Radial component of motion: 


i Ou v Ou ie OP ry 1 OPro il 
o(u e+ 2.55 - 2) = py. Suey 00 +5 Pr, Poo) 
Op Ou 1 du u 02v 1” ag 
Se eA ae he by vt > Ov0e oa 
OA / Ou Uu 1 30 0? Ou Ou 
ia aoe ae 26;) chat Harare ere 
lO ee 0? thay 
so (hear ea pee) 
1 Oe ie row ov v 
+ +08 ae 00 z] 
1 Ou I oh u 
A Y 2u (5 ry 00 “a 


Cross radial component of motion: 


ov v Ov Uv OD,6 1 dfoe 2 
o(u a + “aq 8 ) ore Pro 


y y¥ Ov x 00 Y 
y 1 02u 1 Ou 02u it Ov v 
: H(; Mr ere mr or is 
Ou fl ou ov v I ateye 
+35 00 + 80 7) Fig 


A 02u eee 1 oy 
+ le tae io tees 
ve M (5 ; 02y 2 ou 
nr \7 Oo" ry 00 
1 Tou ee ae 
+ ag (ae 8 at ea 
Le Of a Za” 2 
ae 0) fe ou u) 
1. SOD Ov 
2 ‘ a = 
Bowe: oe) 


(SIE, 
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Continuity equation: 


Energy equation: 


oleae + Foy) — (4 Sy 


poe 4 0 (ee 04 1 0 Le oT 
ae 


Pr ° eal tis 
which reduces to 


e0 (« ahs 4 v oF) (u op fog Sy Sh) 


Y y 00 or ry 00 
a 1 0 a oT iL 0 Lh oT °) 
?. lane ae aes (aS vp hus 
by putting 
t=c,T, 


where ch stands for the dissipation function given by 
Ach =A (¢,, + egg)” + 2 (6, + eg + 2 76) w 
= (A+ 2p) (2, + ef 4) + 2(e,, egg + me) - 


With the usual boundary layer approximations [3], we assume the following mag- 
nitude of the flow parameters within the boundary layer: 

(a) v is of the order (1); 

(b) 6 thickness of the boundary layer is very small compared with 7. Because 
of (a) and (b) it follows from the equations (8) and (7) that uw is of the order of 6 
and yu (as well as 4) is of the order of @ 6%; consequently » is of the order of 6? or 6 
of the order of »1/?. Then equation (8) reduces to 


0 1 0 
A 2 a irae 7A A ae . (10) 
Equation (6) reduces to 
eae Dye (11) 
oP ore 


i. e., the variation of pressure p within the boundary layer is of the order of @ 6 and 
can therefore be neglected. Hence p inside the boundary layer will be given by the 
pressure outside the boundary layer, namely by the pressure in the Prandtl-Meyer 
expansion in the non-viscous case py) = fo(8). Equation (7) reduces to 


ov BOO! ea =A dv 12 
e(u +e selon do. or bee te) 


Dissipation function ch becomes p(dv/dr)*? and equation (9) reduces to 


Or i Gee eS | ae es) 13 
{us +2 se) %—- > ae rw lg am Bren fee Ue) 
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The viscosity-temperature relation [4] is written as 


where 


The equation of the state of the gas is 


pia Onl: 


.ZAMP 


(14) 


(16) 


Thus we have equations (10), (12), (13), (14) and (16) for the determination of the 
five unknown quantities namely, u,v, T, @ and yw. The pressure p is given by its — 


value in the non-viscous case. 


Representation of the Temperature Profile 


We suppose that T = T(v) is a solution [5] of equation (13). Then equation (13) 


transforms into 


Py or 


Fon (enemas ye 0 (2g 4, a ar by 


Pe ai aah. 7S NAOT” Hr 


by neglecting the variation of y within the boundary layer and observing that 


ee He OE cate Ate Bete 
te jvm” One e 0r 06 to 36° 


Utilizing equation (12), the above equation further reduces to 


1 dpb, 1 0 dv 
< + Fgh Cp Ty + 0) + Ty ty (Pr = I ee (x ¥) 


= Fe 6 Ton + 1) ( 


Py 


Assuming Py = 1, we find that c, T, + v = 0 satisfies the equation, i. e., 


if 
EDU ik aa (v§ — v?) 


p 


gives the variation of temperature within the boundary layer. 


Boundary Layer Momentum Integral 


Integrating equation (12) between the limits y = Vy tor =r, + 6, we have 


Ywrto ony 


Tw 


i ou v ov dpy 
fe(w og caval ea [re]2 Or uF ) ar 


2 
ae 


Eo EEE 
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The first term of the integral from the L.H.S, is 
ryt rw+d 


ov ing, 0 
jew oy ~~ Mavet fv 99 (8%) ar 
Tw 


Tw 
and the second term of the integral from the R. H. S. becomes 


rwr+d 


‘ ma =| ine dv 
| or ‘ Or , ee ge cae 


Tw 


by utilizing the boundary conditions u = v = 0 when = ,, and u = uy 
du : 
U=U, (Se) =0O when v=y,-+ 06 and equation (10). 


Neglecting the variation of y within the boundary layer, the momentum integral 
reduces to 


1 vier ta d Twt+o O 
¥ 1 2) F v 
ia = ae Jeera y2 ae 4 0 
Qo Yo My + Ten / 00 (2 ) dy Vw dé / dy Hol Or ye 
Tw Tw 
We define 
rwtd 
: woe uv 
o** = momentum thickness = / ( 1— ) ar , 
20 Yo ry) 
Tw 
rytd 
6* = displacement thickness = / (: eae ) dy . 
20 Yo 
Tw 


Further from the non-viscous motion outside the boundary layer, we have 


do ie apg 
20 Yo (Sp + %)= — ae. 
Substituting the value of dp,/d0, the momentum integral finally is transformed into 
ao=* 1 “du, Ott 20, Ug Ly ( dv ) 
- -— (2 6** + dF . 6 = - 118 
dé? a. ao aloe bi Qo dO Pa Uo Op UE NOY a 


We apply the Polhausen method to solve the above equation. We introduce the 
variable 7 defined [6] by 


such that 7 = 1 when r =7,, + 6, and 7 = 0 when r= 7“y. We assume v/v, = ®(7) 

where @ is a polynomial in 7, with the boundary conditions v = u = 0 for 7 = 0, 

and v= v,, D = 1, 6” = 0, for 7 = 1. ; 
Assuming at first a parabolic distribution of velocity, 1. ¢., On) = 2a ae 


1 
0) 
bt = 6 — 6 [O(n dn = 5 
0 


ZAMP X/38 
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and ; 
2 
ge = 6 [ (y) (1 ~ Pn) dn= zs 9 
0 
But 
ov Q 
pote |) =, 2 a ; 
( dl Yo ( ) 6 Qo 
and hence 


Therefore, equation (18) becomes 


do? 2 9 dvy 4A 2 : dQ 15 =} Ae 30 My Ow Nw 19) 
de® vA eae J) ae ath ( 


From the solution of the non-viscous [7] case, we find 
Uo = Imax SIDAO, Uy = A Qmax cosh 6, 9 = Q59 COS*A O, (20) 
ox Pso cos"P*? 4 6, Zo a 10 cos? 6 5 iii a sec2F +14 7) ; 


where 


= (2 8 eel 


oa == © for (air) == 145 
and 


: a 
Tan = Mi + Oh = uh + 608. 


Equation (19) with the help of the relations in (20) simplifies to 


do? 
dé? 


, 


+ 8 [— (944) Atando + > tan ni 6] = [30 Fo Hy w 
: 05 Yo 


which is a linear differential equation whose solution is given by 
ee E [23+ [30 “#8 w oe Ow cos(9+48- 15/A?) 16 ao] cos —(9+4B- aad 6, (21) 


63 being the integration constant. 


In the present case of flow with pressure gradient [8], @,, /4,, is not constant but 
is given by 


Hoo P ") 2B +2 
= 1G, pest 
iy Op ( RT cos A6 
{from equations (14) and (16)], where Cisa constant in the case of constant surface 


temperature. Replacing r,, by 7,,,sec??+1, 9 (equation of the free stream line), 
solution of 62 becomes 


C Meo P 
=| 62 30 we + Meo Pas i ; 
+ RTx C50 2° Iman [eos ane au 2) aig 2 
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[after substituting the numerical values of J and fin equation (21)], or 


; Cc 
Oo. = [8 aK Yws © Hoo Pso 
0 130 Riot, Pa... 7(0)} cos*4A 6 (22) 
(after integration), where 
sind 0 
Be) = 68 cos®84 6 
. G2 i 67-65 67-65-63... 3 
ioe ieee 2 pik Einar 4 “A : 
“ | pee 90 ee Ge ne ee eee cost 0] 
67- 65+ 63-.- 


~ 3.1 mu AO 
68-66-64... 4.9 stan (F 0 =): 


The function /(@) is represented graphically in Figure 2. 


Figure 2 


If we now assume the Mach number of the incoming stream to be unity, i. e., 
Ma, = 1, the initial Mach line will be perpendicular to the first wall (Figure 1d) 
and 6, is the thickness of the boundary layer at the point O, on the first wall; and 


its value [9] is 6) = 0-632 |/v,9 x C/c*, x being the distance of the point O, from the 
leading edge. Further, the sheering stress at the point O, is 


Eso c* c*¥ C 
Te 02. | ge 
and is same as that from the right hand side obtained by putting 0 = 0, i.e., to 
the value 


2. C* boo PsoC 
89 Qs0 R Too” 


6=0 


Tw 
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Thus, Peo Hoo! Too = 0°095 Myo Oso, Where the suffix so represents values at the - 
outer edge of the boundary layer. This reduces equation (22) to 


C 
52 = [0-399 Yeo ae + 72 099% use /(0)| cos, 6. (23) 


However, we can rewrite equation (23) as 


/2 
sie [: pies (0), cos’/24 6, (24) 
Oo x 
where 
i = Tes Sws , (25) 
x Hi 


Relation (24) is shown in Figure 3 for different values of y and 6. Further we 
assume that 


Vw 


2 
gZe Cole (GF 5 (26) 


o being the deviation factor by which the Prandtl-Meyer stream line differs from 
a circular arc. Its appropriate values for different Maj (i. e., for different g) can 


1] 


10 


0-9 


08 


07 


0-6 


ie 


Figure 3 


be determined from the equation of the stream line. The limiting value of o in the 
case of a circular arc is unity. 


Taking h = 6,/6 and 


5, = 0632 |/oC— DE 


c* 


and replacing 6 by w, we obtain from (24), (25), and (26), for a given expansion 
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May (i. e., y) the coordinates of the virtual point O as given by 


(=z. _|/ Re (1 — h? cos714 w) 
6,/ 0-632 \ gar Ee 2 
5 [« — h® cos?4A wm) + 17-5 h® f(w) cos™4A w cot” iB 

(27) 


uf 
x i aa ? 1/2 , 
[17-5 h? f(@) cos™4A w cot o| 


ay _ Pp 
(“2*) ba 1 |/ Re (1 — h® cos™1A @) cot o 


[a — h*? cos™4A w) + 17:5 h®? f(w) cos74A w cot a 


2 
(28) 


1 
x ea 


[: 7:5 h? f(w) cos74A @ cote o| 


1/2 ? 


where 1 < / < sec’!’*4 9. The relations (27) and (28) are numerically evaluated for 
different values of Maj. 

From PRANDTL-MEYER’s expansion, it is easily seen that if Ma; = 1-093, then, 
gy = 1:4°, wm = 26°, o = 1-1 [10], f(w) = 0-3 (from Figure 2), h = 1-84 (from Figure 3, 
for y = 2, say); and similarly for other Maj. The typical Reynolds number in the 
case of laminar flow without separation [11] are 10° and 10°. In the following table 
the value of REyNOLDs’ number has been taken to be 10°. 


May 1-03 1-06 1-093 
sale 1-06 2-06 4.5 
do 
“ws 3-3 10% | 4-22%102 | 4-75 x10? 
0 


Complete Solution of the Expansive Flow 
D> Uo (2 fee 7?) = 4 mas cosd 6 (2 L/D 7?) » P=Po= Pso cos"F +? 4 6. 


From equation (16) 


p Peo cos"? tA 0 
= = 2 
ae Ug SS yee | 
R [T+ 2G (1-2-7? ) 


From equation (14) 


Tsp ites v9 1-2-7 )| 


p 


_ From equation (12) 


2U9P 
2 0) 0 1 2) 4 ZP234 Bein 1 6 = Cf Oe. 
“* 209 (l= 1) @ Fe ea Ree Too 6 0 


2 
-++ se Uo (2 U UI” Imate sind 0| ‘. 
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Further Calculation with a Different Velocity Profile 


Here we assume a different velocity profile, namely u/V9 = P(n), where 
@ = Cin + Cyn? + Cz 7%, a third degree polynomial in 7, which becomes 
Din) = 29S Oe ls 


under the boundary conditions 7 = 0, ® = 0 and n 1,@ 1, d’ = 0 and 8 ==nOe 
Starting with this velocity profile, the differential equation (19) for 6, in the 
present case, becomes 


dé? 20.) ele dU, 2 dQ ie 56 He = BOVE Oarilen 
dé? 3 Wy A Gone 3 02 Uy 


? 


the solution of which can be written as 


C iG 
Sites [0-399 M59 ye + 9°563 7(6)| cos?95/2] 6 , 


where 
r sind 0 
19) = | 393 | 
—— cos8?2/3 4 6 
ers ‘ 
319 319+ 313 315. Gi se 
esi 2 Se ee 4 a ST eres So leo 108 
x [1+ 360 so , SIRI Coe ae eae ee A 6| 


310 -3139.507s- Gel 3 
) 322 Sorat 


- [costa ae A d ee 0 a os ‘ ees 6 Ri 00] 
The complete solution of the expansion in the present case is given by 
V = Uy (3 7 — 3 4? + 9) , 
P = Pyy cos*F +2] 6, 
Peo Cos?Ft? A 6 


2 
v9 


- a “<i eke ee eee SD 
R [T+ e— (1-30-39 +H) 
D 


> 


_¢ Ho |r v5 (Hose ae 

Be | To 2C, 39 = 3 7? + 4 )]. 

eet E (2B + 2) Ap,» cos?*+14 9 sina 9 — C Hoo, 9 ¥oPo (1 — n) 
3% (1—7)?@ lre a Teer; 


22 ; 
+ > Y% (3.7 — 3.9? + 7°)? Iman Sind 0| : 
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Zusammenfassung 


Die Wirkung der Reibung auf die Prandtl-Meyer-Expansion wird untersucht 
mit Hilfe eines einfachen mathematischen Modells und mit der iiblichen: Grenz- 
schichtapproximation nach POHLHAUSEN. 


(Received: January 2, 1959.) 


Eine Klasse inhomogener Bessel-Gleichungen 


Von VAcLAv VoptéKa, Pilsen, Tschechoslovakei 


Der als Beispiel gewahlte Typus von Gleichungen enthalt als einen Spezialfall 
auch die gewohnliche Gleichung fiir modifizierte Besselsche Funktionen. Die tech- 
nische Bedeutung unserer Ergebnisse soll am einfachen Problem aus der Warme- 
leitung erlautert werden. 


1. Allgemeines 


Es sei 
ae) dy i 
L,{y) = FPR P(*) pe q(*) ¥ (1) 
und f(x, w) irgendeine Funktion von # und o. 
Setzt man 
gle, 0) = [Lp Lee, odo, (2) 
so ist fiir beliebige reelle #), # der Ausdruck 
o 
F(x; 9, 8) = te f(x, o) do (3) 
Bs 


augenscheinlich ein partikulares Integral der inhomogenen Gleichung 


L(y) = g(*, 8) — g(%, Bo) - (4) 


604 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


Die Anwendung dieser fast trivialen Bemerkung soll nun am Spezialfall 


d*y ik ay 
Ly) = Ge +a y (5) 


gezeigt werden. 


Beispiel 7. Unter Zugrundelegung von (5) setze man 


f(¥, o) =sinvxo@. 
Eine kleine Rechnung gibt dann nach (2) und (3) 
g(x, w) = S [(42 + x2 wm? — 1) coswox—wxsinw *], 


1 


TEC Oy) O)) = i (cos ®) « — cos dx) , 


so dass zum Beispiel die Gleichung 


L(y) = = [(m? 2? + x2 — 1) cosa x — axsina x — #4 1) 
den Ausdruck (I), Ay sind natiirlich modifizierte Besselsche Funktionen, C,, C, 


beliebige Festwerte) 


y = Cy Ig(x) + Cy Ko(¥) + - sin 


zur allgemeinen Lésung hat. 

Von diesem verhaltnismassig kiinstlichen Spezialfalle, welcher nur aus metho- 
dischen Griinden eingefiihrt worden ist, gehen wir nun zu mehr sinnvollen Erwa- 
gungen tiber. 


Beispiel 2. Es sei — natiirlich stets mit (5) — 


Mes o) = po EE® 
und damit 
L, [f(%, w)] = = (cos wm — # sin? w) e*°S? = 2 : Lee (sin w e% 008.) 
% x dw : 
Nach (2) und (3) haben wir in diesem Faile 
1 By 
g(¥, w) = Gy Sing ef °F (x; Bo, 8) = [etc dw, 
Do 
und der Ausdruck 
o 
y= Cols C. Kia) ee de (6) 
Bo 
erscheint damit als allgemeine Lésung von 
eal 
L,{y) = — (sin 6 e*°°8? _ sin 9, e* 08%) | (7) 
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Danach hat man zum Beispiel in 
uv 


y = Cy I(x) + Cy K(x) + / eFC8® Gey 
0 
die allgemeine Lésung der Gleichung 


1 Z 
L,(y) eS = sin @ e* © o , 
Es gibt Werte von #, #, fiir welche sich (7) in die homogene Gleichung der 
modifizierten Besselschen Funktionen verwandelt. Dies gilt offenbar fiir 
C= wi OSes nh ED De aan, 
und die Tatsache 


nm 


ee Bey = (n —m) x Tala lr) (oe Oy Sire eZ on. 


mm 


steht damit im vollen Einklang. 


Beispiel 3. Es sei nun 


und 
f(x, ) = sinh(¥ cosh@ — v @) . 


Es ergibt sich 
s a 
L, (f(#, o)] = (sinh? o _ =) sinh (¥ cosh m — vy w) + s cosh w cosh (¥ cosh w — ¥ w) 


= Db 4 sat [(~ sinh w + v) cosh(4 cosh w — ¥ @)] 
x2 da 


und nach (2), (3) ist daher 
g(% wo) = ae [(x sinh w + v) cosh(¥ coshw — v @)], 
x 


0 
F (x4; 85,8) = | sin coshw — vw) do. 
Oo 


Wir gelangen so zur Erkenntnis, dass der Ausdruck 


0 
y = C, I,(#) + C, Ky(%) + [sinh (w cosh w — v m) dw (9) 
Io 


606 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 
die vollstandige Losung von 

L(y) = = [(x sinh # + v) cosh(x cosh B — » B) (10) 
— (x sinh #, + v) cosh(¥ cosh Jy — » Bo) ] 


darstellt; L, erklart sich nach (8). 
Im Gegenteil zu (7) enthalt der Typus (10) nicht die zugehérige homogene 


Gleichung L,(y) = 0. 
Unser letztes Beispiel gehért zur Theorie gewohnlicher (das heisst nicht modi- — 
fizierter) Bessel-Gleichungen. 


Beispiel 4. Es sei 


ay et ad 2 
Ld) = perch oat beeen oy 


und wahlen wir in naheliegender Weise 

f(x, w) = cos(#’sinw — vw). 
Damit ergibt sich 

2 


v : Te a5 : 
Les ON = (cos? o— =) cos (¥ sinw — » w) — =| sine sin (+ sinw — y @) 


d 2 : 
Sp OVER [(¥ coswm + v) sin(¥ sinw —» o)], 


und wir erhalten 


if 
g(x, w) = ox [(¥ cosw + v) sin(¥ snw —yo)], 


o) 
F(x; Bo, ) = | cos(x sino — yo) da. 
Bo 


Es ist daher 
o 
¥ = Cy Jy(4) + Cy Yo) + [cos (v sinw — » w) dw (12) 
i 
die allgemeine Lésung der Gleichung 


1 
L(y) = =e [(¥ cos # + ») sin (x sin ® — » ) (13) 
3 
— (¥ cos&> + ») sin(x sind, — v J,)]. 


, Der Leser gebe sich bitte Rechenschaft iiber die zahlreichen Spezialwerte von 
ee y ae ia sich (13) in oA homogene Bessel-Gleichung L,(y) = 0 verwandelt. 
le Konsequenzen fiir die zugeh6ri i i 
mise acon: gehorigen, in (12) vorkommenden partiku- 
Unsere Betrachtungen liessen sich offenb 1 

: : ar zu einem System ausbauen und 
auf ore (nicht notwendig Besselsche) inhomogene Gleichungen erweitern, wobei 
man aut ganz neue, merkwiirdige Tatsachen kommen wiirde. Jedoch ist das an 
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dieser Stelle nicht beabsichtigt — das Ziel dieser Aibeit ist es ndamlich, nur eine 
Grundlage zur mathematischen Behandlung gewisser vom Verfasser verfolgter 
Fragen aus der Warmeleitung zu schaffen. Indessen sei im folgenden Abschnitt 


mindestens eine elementare Frage aufgeworfen, wo unsere Ausfiihrungen zur 
Geltung kommen. 


2. Momentane Warmequelle langs eines Kreisbogens 


Es spiele sich zur Zeit ¢) langs des Kreisbogens 


€—x7 = Reosyp, n—y =Rsiny, €—-29=0, ag SpSe + (14) 


’ 


eine momentane Warmeexplosion von konstanter Starke o ab. In einem allseitig 
unendlichen, homogenen und isotropen Medium mit physikalischen Festwerten A, c, 
oe entsteht dadurch das Temperaturfeld 


& 


Ro h(t; to) — Jets te—el0l)8s R/4 gt Ge ft een 
ul, ¥. 234) = gos osm Gain [w* +( We es ed ry V0) dy: 
&o 
A at Rw 
@=-—, A(t;t)=-— > [1+sgn(t¢—)], «= =, 
tp Mise abt 2 a? (t — t) a 
= 2 4 tek ee (ee OP (0))2 
Baie Sy ee at ae) 


Durch Einfiihrung der neuen Integrationsvariablen w = y'— po, dw = dy er- 
scheint unsere Temperaturfunktion (15) in der Form 
E— Wo 


u(x, y, 23 t) = Ro hit; to) gle sai")? Ree eX COSO Toy (16) 


8 ae 2 a xa ty) 3/2 
€y9— Vo 
Ist y, konstant, das heisst, interessiert man sich um die Temperaturverhaltnisse 
in irgendeiner Ebene des Biischels mit der Achse x = +), y = y), so stellt nach 
(6) und (7) das bestimmte Integral aus (16) eine partikulare Lésung der inhomo- 
genen Bessel-Gleichung 


d*y 1 dy 
dx? x dx 


y = — [sin (yy — 6) e208 sin (6 = ye) OEM], (17) 
M4 


dar. So kommen unsere obigen Erwagungen zur Geltung. ar HM 
Im allgemeinen ist jedoch das Temperaturfeld (16) ungewohnlich kompliziert. 

Nur im Falle ¢ = ¢, + 2” kommen wir auf die elementare Formel?) fiir eine ring- 

formige Warmequelle, wahrend sich schon eine Halbkreisquelle mit ihrer Tempera- 

turfunktion 

Ro hit; ty) [w+ (2210)? + RP]/4.a (tty) 


u(x, ¥, 2; t) = 8a? m2” (t — t,)32 


eee (18) 
SCuaeel, Gee 2 | sinh (x cos w) dw 
F ‘ 
als recht kompliziert erweist. 


1) H. S. Carsvaw und J. C. JAEGER, Conduction of Heat in Solids (Oxford 1947), S. 220. 
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Der Leser mége selbst das Studium der entsprechenden Gleichung (17) 
vornehmen und sich namentlich mit dem Temperaturzustand in den Ebenen 
Wo — 9 = 0, 2/2, m naher befassen. 


Summary 


New particular solutions are given for certain types of nonhomogeneous Bessel 
equations. Previous considerations allow of considerable extension and their 
technical importance is illustrated by an example in heat conduction. 


(Eingegangen: 4. April 1959.) 


Varia — Miscellaneous — Divers 


6th National Symposium on Reliability and Quality Control 
in Electronics, Statler Hilton Hotel, Washington, D.C., USA. 


In Washington findet vom 11. bis 13. Januar 1960 das 6. National Sym- 
posium on Reliability and Quality Control statt. Der Zweck dieser Tagung, welche 
vom Institute of Radio Engineers (IRE) und anderen wissenschaftlichen Ver- 
einigungen durchgefiihrt wird, ist, Spezialisten, welche auf dem Gebiet der Be- 
triebssicherheit und Qualitatsitiberwachung iiber Erfahrung verfiigen, zusammen- 
zufiihren. In Vortragen und Diskussionen wird ein ausgiebiger Gedankenaustausch 
gepflegt. Programme und Anmeldeformulare sind zu beziehen bei C. M. BEYER, 
OASD, Room 3D1031, The Pentagon, Washington 25, D.C., USA. 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Elementare Informationstheorie. Von Hrinz ZEMANEK (R. Oldenbourg, 
Wien und Miinchen 1959), 120 S., 28. Abb., 2 Tabellen. 

_ Der Verfasser versucht im vorliegenden Biichlein die Informationstheorie in 
einer dem Ingenieur zuganglicheren Form darzustellen. Er unterscheidet zunachst 
2 Hauptprinzipien, welche getrennt betrachtet werden kénnen. Zunachst wird das 
logische Netz entwickelt, das heisst die Form, in welcher sich «Information» aus- 
driicken lasst. Dazu gehéren die verschiedenen Kodierungsméglichkeiten. Eine um- 
fangreichere Betrachtung ist der Information als statistischer Vorgang gewidmet. 
Die Entwicklung des optimalen realisierbaren Kodes, zum Beispiel mit Binar- 
Zeichen, wird auf geometrischem Wege in Abweichung zu SHANNON entwickelt, 
wobei die angegebene Methode eine kleinere Redundanz ergibt. Beispiele iiber In- 
formationsgehalt werden aus der Sprache und der Biologie besprochen. Ein drittes 
Kapitel ist der statistischen Zuordnung gewidmet, welche eine Rolle spielt bei der 
gestérten Ubermittlung, das heisst bei der Berechnung der maximalen Kanalkapa- 
zitaét eines Ubertragungsvorganges bei Anwesenheit von Stérungen. Der Ubergang 
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auf kontinuierliche Vorgange und eine kurze Diskussion der Quantisierung eines 
kontinuierlichen Vorganges beschliessen den Text. Angefiigt ist ein recht vollstan- 
diges Schrifttumsverzeichnis und 2 Tabellen. 

Das Biichlein bietet auf engem Raum eine sehr gute Einfithrung in die Informa- 
tionstheorie. Etwas stérend scheinen zunachst die zahlreichen neu eingefiihrten 
Begriffe mit ihren zugehGrigen Zeichen zu sein, sie erleichtern aber doch begrifflich 
die Trennung des Informationsgehaltes vom Entscheidungsgehalt eines gewahlten 
Kodes fiir eine bestimmte Informationsquelle. H. WEBER 


; Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLUGGE und E. Creutz. Band 45, 2. Teil: Instrumentelle Hilfsmittel der Kern- 
Physik (Springer Verlag, Berlin 1958). 544 S., 293 Abb.; DM 128.—. 

In dem in englischer Sprache abgefassten Band vermitteln 12 wohlbekannte 
Autoren ihre Erfahrungen auf dem Gebiet der Detektoren und Nachweismittel fiir 
geladene und ungeladene Teilchen sowie fiir elektromagnetische Strahlung. Die 
Wechselwirkungen von Partikeln und Strahlung mit Materie werden an verschie- 
denen Stellen des Buches diskutiert. Dem Gebiet der Gaszdhley sind drei Kapitel 
gewidmet, die sich mit Ionisationskammern, Proportional- und Auslésezahlern 
befassen. Nach einem theoretischen Uberblick iiber die Vorgange im Kammergas 
werden eine Fiille von Einzelheiten iiber die technische Ausfiihrung solcher Mess- 
gerdte mitgeteilt. Spezielle Anforderungen an Detektoren dieser Art, wie sie zum 
Beispiel zum Nachweis sehr schwacher Intensitaéten, zum Zahlen von Gasproben, 
Messung von Teilchenspektren usw. gestellt werden, sind ausfiihrlich diskutiert. 

Das Kapitel iiber Szintillations- und Cerenkow-Zahler ist relativ knapp gehalten. 
Ein Uberblick vermittelt Kenntnisse iiber die Wirkungsweise und die verschiedenen 
Konstruktionen von Photomultipliern mit Beriicksichtigung der zurzeit handels- 
ublichen Typen. Die wichtigsten Kristalle, wie NaI, Lil, CsHal, sowie einige orga- 
nische Phosphore werden im Hinblick auf ihre Herstellung, Lichtausbeuten, Zer- 
fallszeiten usw. besprochen. Edelgase als Szintillatoren werden dabei nicht beriick- 
sichtigt. Auch vermisst man einen Hinweis auf diejenigen Vorgange, die fiir die 
Konversion der im Kristall dissipierten Energie in Licht verantwortlich sind. Tech- 
nische Probleme, wie sie sich bei der Kombination eines Phosphors mit einem 
Elektronenvervielfacher stellen, so zum Beispiel die Ausbildung von Reflektoren 
und Lichtleitern, werden kurz berithrt. Bei der Behandlung der Anwendungs- 
méglichkeiten sind einige neuere Methoden, wie die Kombination eines schnellen 
und eines langsamen Phosphors zu einem sogenannten Phoswich oder die Ausniut-’ 
zung der Abhangigkeit der Zerfallskonstanten einiger Szintillatoren von der 
Tonisationsdichte, nicht erwahnt. Der Abschnitt iiber Cerenkow-Zahler wird durch 
eine alles Wesentliche enthaltende Charakterisierung der Eigenschaften der Ceren- 
kow-Strahlung eingeleitet. Bei der Diskussion der spezifischen Eignung des Ceren- 
kow-Zahlers zum Beispiel als Schwellendetektor, Zahler zur Bestimmung von 
Teilchengeschwindigkeiten, Teilchenladungen und Flugrichtungen wird dem 
Gebiet der Hochenergiephysik besondere Aufmerksamkeit gewidmet. Im Kapitel 
iiber den Nachweis von Neutronen werden diejenigen Streuvorgange und Kern- 
reaktionen besprochen, die zur Registrierung von Neutronen benutzt werden. Die 
apparativen Hilfsmittel, wie Gas- und Szintillationszahler, Wilson- und Blasen- 
kammern, photographische Emulsionen sowie die Aktivierungsmethode, werden 
hinsichtlich ihrer spezifischen Brauchbarkeit bei diversen physikalischen Problemen 
untersucht. Aus der Fiille der Anwendungsmethoden werden die Messung von 
Neutronenfliissen, die Bestimmung von Quellenstarken, die Neutronendosimetrie 
und die Aufnahme von Energieverteilungen besonders hervorgehoben. Auch an 
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dieser Stelle werden Probleme der Hochenergiephysik separat behandelt. Im Hin- © 
blick auf die Wichtigkeit der Koinzidenztechnik in bezug auf sémtliche mit Detek- 
toren zusammenhangenden Probleme wird dieser Methodik ein separates Kapitel 
gewidmet, in dem zahlstatistische und elektronische Fragen erlautert werden. 
Theoretische Betrachtungen iiber die durch ionisierende Teilchen in tibersattig-: 
ten Dampfen ausgeléste Trépfchenbildung und uber die Erzeugung einer Uber- 
sittigung durch adiabatische Expansion leiten das Kapitel tiber Nebelkammern ein. 
Anhand der Beschreibung einiger typischer Expansionskammern wird der Leser mit 
zahlreichen technischen Einzelheiten vertraut gemacht. Zur Kenntnis der Diffu- 
sionsnebelkammer vermitteln die Autoren einen Abriss der Theorie dieses Instru- 
mentes und die Beschreibung der praktischen Ausfiihrung einer solchen Kammer. 
Die Aufzihlung der zahlreichen und wohlbekannten Anwendungen ist stets von 
wertvollen Hinweisen auf die Auswertung und die Interpretation der Messresultate 
begleitet. ; 
Die Blasenkammer beschreibt der Erfinder dieses wertvollen Instrumentes in 
einem kurz abgefassten und doch sehr reichhaltigen Kapitel, in dem er sich mit den ~ 
Vorgangen in der Kammerfliissigkeit und den verschiedenen Konstruktionsméglich- 
keiten auseinandersetzt. Wie an anderen Stellen des Buches fehlt es auch hier nicht 
an Hinweisen auf zahlreiche Anwendungsméglichkeiten. : 
Ein recht umfangreiches Kapitel beschaftigt sich mit der Anwendung der Photo- 
emulsionstechnik in der Kern- und Partikelphysik. Die Vorbereitung der Emulsionen 
zur Exposition, die Entwicklungsmethodik, das mikroskopische Auswertungs-. 
verfahren und die physikalische Interpretation vonSpuren werden detailliert be- 
handelt. Daneben findet man eine Fiille von numerischen Informationen, wie etwa 
Energie-Reichweite-Beziehungen, Schrumpfungsfaktoren usw. : 
Die aussere Gestaltung des Buches entspricht dem Habitus der bekannten — 
Fliiggeschen Handbiicher. Es umfasst 544 Seiten und enthalt 293 Figuren. i 
Es ist erfreulich, dass auch der Experimentierkunst des Physikers im Rahmen 
dieser Handbuchreihe gebiihrend Rechnung getragen wird. F. HEINRICH. 


Physikalische Dynamik der Atmosphare. Von ALBERT und FRIEDRICH 
DrFant (Akademische Verlagsgesellschaft mbH., Frankfurt a. M. 1958). 527 S., 
139 Abb.; DM 60.—. 

Unter den wenigen deutschen Lehrbiichern der theoretischen Meteorologie 
fehlte bisher ein Werk, das planetarische Probleme, wie die allgemeine Zirkulation ~ 
der Atmosphare und die Grundlagen der numerischen Prognose (Vorausberechnung 
des Str6mungsfeldes mit Hilfe elektronischer Rechenautomaten), behandelt. Diese 
Liicke schliesst das vorliegende Buch der beiden dsterreichischen, als Forscher und 
Hochschullehrer auf dem Gebiete der Meteorologie und Ozeanographie bekannten 
Autoren Prof. A. und Fr. Deranr (Vater und Sohn). 

é Der Titel ist offenbar in Anlehnung an die Physikalische Hydrodynamik von 
Vi BJERKNES gewahlt worden und weist auf die innige Verflechtung verschiedener 
PE ve easelies Disziplinen bei der theoretischen Deutung atmospharischer Vorgange 

in. 

_ Der Stoff ist itbersichtlich gegliedert. Bis zum 7. Kapitel entspricht der Aufbau 
cinigermassen dem Herkémmlichen: Statik, Thermodynamik, Stabilitat vertikaler 
Luftsaulen, Strahlung, Kinematik, Bewegungsgleichungen auf der rotierenden 
Erde, einfache Horizontalstrémungen. In den folgenden Kapiteln werden neuere 
Begriffe behandelt: Zirkulation, Vorticity, Diskontinuitatsflachen, Turbulenz, Aus- 
tausch und Reibung. Lebhaft zu begriissen ist sodann die eingehende Darstellung 
der atmospharischen Energetik, der allgemeinen Zirkulation und der numerischen 
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Vorhersage. Diese neuen, in Entwicklung begriffenen Teilgebiete lassen sich zurzeit 
noch kaum in ein geschlossenes Lehrgebaude der atmospharischen Dynamik ein- 
ordnen ; doch haben sich die Autoren bemiiht, einen umfassenden Uberblick iiber den 
heutigen Stand der Theorien zu geben. Das letzte Kapitel ist den atmospharischen 
Gezeiten gewidmet. 

Der Auswahl des Dargebotenen ‘darf man — von Einzelheiten abgesehen — be- 
friedigt zustimmen. Die Formelsprache bleibt durchwegs elementar. Leider haben 
sich in Text und Formeln einige Vorzeichen- und andere Fehler eingeschlichen, die 
bei einer Neuauflage ausgemerzt werden sollten. Der sprachlichen Formulierung ist 
unseres Erachtens zu wenig Sorgfalt geschenkt worden; einige Satze wirken geradezu 
grotesk. 

Trotz dieser Schénheitsfehler ist dem Lehrbuch wegen seines umfassenden In- 
haltes eine weite Verbreitung an meteorologischen und geophysikalischen Insti- 
tuten, bei Fachleuten und Studierenden zu wiinschen. Der Verlag hat fiir einen sehr 
gediegenen Druck gesorgt. W. KuHN 


Ballistik. Von R. E. Kurrerer, 3. Aufl. (Friedr. Vieweg & Sohn, Braun- 
schweig 1959). 304 S., 176 Fig., 21 Tabellen; DM 38.-. , 

Seit dem Erscheinen der ersten Auflage des Buches von KuTTERER im Jahre 
1942 hat die Ballistik dank der Entwicklung der Raketentechnik und der Elek- 
tronik bedeutende Fortschritte gemacht, und zwar sowohl im Hinblick auf ihre 
Aufgaben als auch in bezug auf ihre theoretischen und experimentellen Methoden. 
Aus diesem Grunde war eine Neuauflage dringend erforderlich und erwiinscht. Der 
Verfasser hat in der kiirzlich erschienenen dritten Auflage den zahlreichen Neue- 
rungen Rechnung getragen und nicht nur die einzelnen Kapitel erweitert, sondern 
auch eine Reihe neuer Kapitel hinzugefiigt, so dass nun wieder ein modernes ein- 
fiihrendes Werk iiber Ballistik vorliegt. Ein erster Abschnitt befasst sich mit den 
Geschossbahnen im luftleeren und lufterfiillten Raum sowie der Pendelung und dem 
damit zusammenhangenden Problem. der Stabilitat eines Geschosses. Sehr ein- 
gehend sind dann die messtechnischen Probleme (Geschwindigkeit, Flugzeit, Rota- 
tion, Flugbahnelemente usw.) behandelt. Die Berechnung und Messung der inner- 
ballistischen Gréssen und die Frage der Leistungssteigerung einer Waffe kommen 
im zweiten Abschnitt zur Sprache: Der letzte Abschnitt-schliesslich geht auf die 
Innenballistik und die einfachern Falle der Aussenballistik von Raketen ein. Das 
Werk von KuTTERER — das mit nahezu 200 Literaturangaben ein eingehenderes 
Studium mancher Sonderfragen erméglicht — hebt vor allem die experimentelle , 
Ballistik und die messtechnischen Probleme hervor. Nicht naher behandelt sind 
dagegen beispielsweise die Flugbahnstérungen, Extremalprobleme bei Raketen- 
flugbahnen und die Ballistik interkontinentaler Raketen. Das groBziigig ausgestat- 
tete Buch wird nicht allein den Ballistiker interessieren, sondern diirfte auch wegen 
den zahlreichen speziellen Messmethoden zur Untersuchung kurzzeitiger oder unter 
extremen Bedingungen vor sich gehender Vorgange fiir manchen Physiker und 
Ingenieur von einiger Bedeutung sein. E. RotH-DESMEULES 


Dynamische Vorgange in linearen Systemen der Nachrichten- und 
Regelungstechnik. Von H. Kaurmann (R. Oldenbourg Munchen 1959). 208 S., 
‘96 Abb.; DM 26.50. . 

Der Autor stellt auf engstem Raum alle modernen grundlegenden Betrachtungen 
an, wie Gewichtsfunktion und Frequenzgang als Kennzeichnung der Netzwerkeigen- 
schaften im Zeit- und Frequenzbereich, Ubertragungsfunktion und Berechnung 
der zugehérigen Gewichtsfunktion, Minimalphasennetzwerk mit den Zusammen- 
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hangen zwischen Phase und Betrag, Bode-Diagramm des Frequenzganges, Stabili-_ 
tAtskriterien, Approximationsverfahren fiir die Analyse und Synthese von Netz- 
werken, Kennwerte von Ubergangsvorgangen, Systeme mit Abtastung. Es wird ~ 
dabei die richtige Tendenz befolgt, das zeitliche Verhalten mit dem F requenzgang : 
parallel darzustellen. Fiir den Ingenieur besonders wertvoll ist die eingehende Be- z 
handlung von wirkungsvollen Naherungsverfahren zur Synthese von Netzwerken, 
sei es auf Grund des zeitlichen Verhaltens oder des Frequenzganges. Erfreulicher- 
weise wurde die besonders in Amerika gepflegte Methode der Polortskurven fiir den 
geschlossenen Kreis eines Regelsystems auf Grund der Pole und Nullstellen des 
offenen Kreises eingehend ‘behandelt. Ebenfalls neu ist das Kapitel Systeme mit 
Abtastung, welche mittels der Z-Transformation untersucht werden. Die Anwen- 
dung eines Digitalrechners im Regelkreis wird dabei kurz gestreift. Tabellen und 
die allereinfachsten Zuordnungen der Laplace-Transformation und der Z-Trans- 
formation bilden den Schluss. Offenbar war beim Abschluss des Manuskriptes das 
Buch von E. I. Jury, Sampled Data Control Systems, noch nicht zuganglich und ; 
fehlt deshalb im Literaturverzeichnis. 
- Die Darstellung des Stoffes ist sehr konzentriert, obgleich die meisten Betrach- as 
tungen anhand von Beispielen veranschaulicht werden. Einem Leser, der sich bereits 
etwas in das Gebiet hineingearbeitet hat, gereicht dies gerade zum Vorteil. Das 
Buch wird daher allen Fachleuten, welche gelegentlich mit diesen Dingen in Bertth- 
rung kommen, ausserordentlich von Nutzen sein. H. WEBER 
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Physics of Meteor Flight in the Atmosphere. Von Ernst J. Opix. Inter- 
science Tracts on Physics and Astronomy, Nr. 6 (Interscience Publishers, Inc., 
New York 1958). 174 S., 14 Abb.; $1.95. 

Das Studium der physikalischen Vorgange beim Eintritt eines Korpers hoher — 
Geschwindigkeit — wie im Falle der Meteore — in die dichtere Erdatmosphare, dessen 
Hauptschwierigkeiten einerseits in den ungewohnlichen Bedingungen und in der ~ 
Tatsache, dass die Beschaffenheit dieser Kérper weitgehend unbekannt ist und 
andrerseits — wenigstens bis vor kurzem — in der Unméglichkeit des direkten Ex- 
perimentes liegen, begann im wesentlichen mit den Arbeiten von Oprx seit 1922. 
Trotz dieser Schwierigkeiten ist es in den letzten Jahrzehnten gelungen, anhand 
verschiedener Modelle dem Phanomen naherzukommen und zwischen Massenab- 
nahme, Geschwindigkeit, Helligkeit eines Meteors, die durch dieses hervorgerufene 
Ionisation Beziehungen abzuleiten, die mit den Beobachtungen verglichen werden 
k6énnen. Im vorliegenden kleinen, aber inhaltsreichen Bandchen von Oprx sind die 
wichtigsten Gesichtspunkte zusammengefasst, und ein grosses numerisches Mate- 
rial fur die verschiedenen Falle ist in den iiber 50 Tabellen enthalten. Die Darstel- 
lung zeigt, dass heute noch keineswegs eine geschlossene Theorie dieser ausseror- 
dentlich komplexen Vorgiénge méglich ist und man vielfach mit der richtigen 
Grossenordnung der Ergebnisse zufrieden sein muss, so dass der Forschung noch 
manche Aufgabe verbleibt. - Die verschiedenen Kapiteliiberschriften lauten: The 
Atmosphere. Classification and Physico-Chemical Properties of Meteors. Meteoroid 
Energy Transfer. Ablation. Details of Energy Transfer and Types of Ablation. Atomic 
Collisions. Meteor Radiation. Some Applications. Die Monographie von Oprk — die 
eine empfindliche Liicke in der astrophysikalischen Literatur ausfiillt — diurfte 
ibrigens nicht nur fiir den Astrophysiker von Interesse sein, sondern auch fiir den 
Ballistiker und den Aerodynamiker, insofern das Problem des Eintrittes eines 
Korpers in die Atmosphare sowohl im Falle der interkontinentalen Raketen als 
auch bei der Wiedergewinnung von Raumraketen von grosser Bedeutung ist. 
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Mitteilungen aus dem Institut fiir angewandte Mathematik — 
an der Eidgenéssischen Technischen Hochschule Ziirich 


Herausgegeben von Prof. Dr. E. Stiefel 


: Nr. 1 


Entwurf eines elektronischen Rechengerates unter besonderer Beriick- 
sichtigung der Erfordernis eines minimalen Materialaufwandes bei gegebener 
mathematischer Leistungsfahigkeit. Von Ambros P. Speiser. (2. Auflage 
1954.) 67 Seiten. Broschiert Fr. 7— (DM 7.-). 

Vorstudien fiir die Entwicklung eines Rechenautomaten mittlerer Grime; 
angepasst an die Bediirfnisse eines mathematischen Instituts. Grundlage 
fiir die elektronische Rechenmaschine ERMETH der Eidgendssischen Tech- 
nischen Hochschule Ziirich. Der Entwurf sieht tooo Elektronenrdéhren und 
300 elektromagnetische Relais vor; als Speicher dient eine magnetische 
Trommel. ~ ae; ‘5 


Nr. 2 


Programmégesteuerte digitale Rechengerate (elektronische Rechen- 
maschinen). Von Heinz Rutishauser, Ambros Speisey und Eduard Stiefel. , 
(1951.) 102 Seiten. Broschiert Fr. 9— (DM 9). Vergriffen. 


. Inhalt: Grundlagen und wissenschaftliche Bedeutung — Organisation und 


Arbeitsweise — Arithmetische Prinzipien — Vorbereitung von Bechegelaney sy 
— Physikalische Grundlagen. 


ae Nr. 3 
Automatische Rechenplanfertigung bei programmégesteuerten Re- 
chenmaschinen. Von Heinz Rutishauser. Nachdruck 1956. 45 Seiten. 


Broschiert Fr. 6.— (DM 6.-). 
Vorstudien fiir die Automation des Programmierens fiir elektronische 


Rechenautomaten. 
Nr. 4. 


An Oscillation Theorem for Algebraic Eigenvalue Problems hs its: 
Applications. By Frank W. Sinden. (1954.) 57 pages. Fr. 7.— (DM 7.-). 


In der Theorie der Eigenwertprobleme bei gewodhnlichen Differential- « 


gleichungen (z. B. kritische Drehzahlen) spielen die sogenannten Oszillations- 
theoreme eine grosse Rolle, welche Aufschluss geben iiber die ‘Anzahl der 
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Knoten der Obarschiyingunisem Der Verfasser aibertrigt diese’ Theorie! auf 


die Differenzenrechnung, womit Oszillations-T. heoreme auch in der numeri-. ‘ 


schen Mathematik ihre Anwendung finden k6nnen. Mg 
ae 
Nr. 5 : . 
Tabellen zur Erzeugung von Funktionen einer und zweier Variablen 
mit linearen Potentiometern. Von Moheb Aziz Abdel-Messth. (1954-) 
33 Seiten. Broschiert Fr. 5.— (DM 5.-). 


Katalog fiir die Erzeugung von rationalen Funktionen durch lineare Poten- — 


tiometer. Es werden alle Schaltungen mit einem oder zwei Potentiometern 


zusammengefasst. Die Tabellen bilden ein Hilfsmittel fiir die Darstellung — 


empirischer Funktionen durch Potentiometer-Schaltungen. 


Nr. 6 
Die mathematischen Grundlagen fiir die Organisation der elektro- 


nischen Rechenmaschine der Eidgenéssischen Technischen Hoch- © 


schule. Von John Robert Stock. (1956.) 76 Seiten. Broschiert Fr. 8.— (DM 8), 
In den Jahren 1952-56 wurde im Institut fiir angewandte Mathematik der 
ETH, in Zusammenarbeit mit verschiedenen Firmen, eine elektronische 
Rechenmaschine konstruiert, deren Grundkonzeption sich in verschiedener 
Hinsicht von anderen Rechenautomaten unterscheidet. Die vorliegende Ab- 
handlung beschreibt die fiir den Bau der ERMETH massgebenden Prinzipien 
sowie den allgemeinen Aufbau und schliesslich auch die Details der Ablaufe 
der Rechenoperationen. 


Nr. 7 


Der Quotienten-Differenzen-Algorithmus. Von Heinz Hiseishates 
(1957.) 74 Seiten. Broschiert Fr. 9.— (DM 9.-). 

Zusammenfassung von drei in der Zeitschrift fiir angewandte Mathematik 
und Physik (ZAMP) erschienenen Arbeiten des Verfassers, nebst einigen 
Erweiterungen und Erganzungen.. 
SSS 
Neuerscheinung: Nr. 8 


Refined Iterative Methods for Computation of the Solution and 
the Eigenvalues of Self-Adjoint Boundary Value Problems. By 


M. Engeli, Th. Ginsburg, H. Rutishausey and E. Stiefel. (1959.) 107 pages. — 


Fr. 17.— (DM 17.-). 
Darstellung einiger bewahrter Verfahren und numerischer Experimente tiber 
die iterative Lésung von Randwertaufgaben bei elliptischen partiellen Dif- 


ferentialgleichungen unter besonderer Beriicksichtigung der Randwert- 
probleme der Elastizitatstheorie. 
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Die Reihe wird fortgesetzt 
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